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Abstract 


We calculate the orbits of the low-energy solar cosmic ray particles in the 
earth’s dipole magnetic field and then estimate on their asymptotic velocity vectors 
before entering this magnetic field. Although the incident region of these particles 
is restricted in the area of higher geomagnetic latitudes than 60°, it is clear 
that there are anisotropic effects, i.e. impact zone effects for the incidence of 
these particles. Furthermore, the precipitation lines of the incidence will be shown 
to be circular around the geomagnetic north pole. 

Next, it will be shown from the comparison of the calculated results with the 
observed data that, in the initial stage of the incidence, these low-energy cosmic 
ray particles tend to impinge upon the earth anisotropically, in accordance with the 
expectation from the calculated results, and that, slightly before the onset of geo- 
magnetic Sc storms, their anisotropic incidence changes into isotropic one, although 
there are some exceptional cases. Finally, criticisms and a possible explanation 
on the features of their incidence after the onset of geomagnetic Sc storms will 


be proposed. 


1. Introduction 


The observations and the detailed studies on the low-energy cosmic rays associa- 
ted with solar flares are being carried on very actively at present (Anderson, 1958; 
Anderson, Arnoldy, Hoffman, Peterson and Winckler, 1959 ; Freier, Ney and Winckler, 
1959; Hakura and Goh, 1959; Leinbach and Reid, 1959; Ney, Winckler and Freier, 
1959; Thompson and Maxwell, 1960). In association with these investigations, the 
studies on the low energy end of the galactic cosmic rays have started, and now the 
knowledge of these cosmic rays is being accumulated (Simpson, 1960). 

The investigations on these low-energy solar cosmic rays made clear in detail 
the relation of solar phenomena to the terrestrial magnetic disturbances, the acceleration 
mechanism of these particles, their propagation into interplanetary space and so on 
(Reid and Leinbach, 1959 ; Obayashi and Hakura, 1960; Winckler, 1960). The incidence 
of the low-energy solar cosmic ray particles to the earth was firstly made clear by 
Bailey (1957; 1959) from the study of the abnormal ionization of the ionosphere as- 
sociated with the cosmic ray intensity unusual increase at February 23, 1956. Then 
during the I.G.Y., many high-altitude observations of these low-energy solar cosmic 
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rays were performed and some properties of these cosmic rays have been clarified. 

In this paper, in section 2., we shall at first calculate in details the motion of 
these low-energy solar cosmic ray particles in the earth’s dipole magnetic field and 
then show these results in some figures. Although a part of this calculation was 
already performed in a previous paper (Sakurai, 1960), we shall investigate here the 
orbits of these cosmic ray particles in more details. In section 3, we shall consider 
the mechanism of the incidence of these particles and feature of the variation of 
geomagnetic fields during the geomagnetic Sc storms. 


2. Motion of the Solar Cosmic Ray Particles in the Earth’s Dipole Magnetic Field 


It seems that the low-energy solar cosmic ray particles, which are ejected from 
the sun associated with the flares and then accelerated in its outer atmosphere, will 
approach the earth while interacting with the plasma clouds or the electromagnetic 
fields supposed to be in the interplanetary space. We may suppose that, before these 
cosmic ray particles arrive near the earth, their orbits will deviate trom the solar 
direction on account of the action of the interplanetary magnetic field because of 
their low energies (Gold, 1959; 1960; Hakura and Goh, 1959; Reid and Leinbach, 1959; 
Obayashi and Hakura, 1960; Sakurai, 1960). According to the studies by Hakura, 
Takenoshita and Otsuki (1958) and Obayashi and Hakura (1960), however, the incidence 
of these low-energy solar cosmic ray particles seems to show such characteristics as 
follows; when these particles impinge upon the polar latitude region of the earth 
before the onset of the geomagnetic Sc storms, there are some anisotropies for their 
incidence, which seem to show that the motion of these particles well agrees with the 
expectation from the theory of geomagnetic effect of cosmic rays (Sakurai, 1960). 
Since some time before the onset of the geomagnetic Sc storms, they will begin to 
impinge upon isotropically, and then come nearly uniformly into, the polar latitude 
region. This change of anisotropic into isotropic incidence appears to mean that, 
according to the present knowledge of the solar cosmic rays, the distribution of the 
low-energy solar cosmic ray particles in the interplanetary space becomes isotropic 
(Meyer, Parker and Simpson, 1956; Obayashi and Hakura, 1960). So, the incidence of 
these low-energy solar cosmic ray particles to the earth shows anisotropy at the initial 
stage and then gradually becomes isotropic, and before long they uniformly arrive at 
the polar latitude region. Hence, we can conclude that, in the initial stage of their 
incidence, these low-energy solar cosmic ray particles will arrive near the earth with 
some anisotropies. From such a point of view, we calculated in a previous paper 
ese 1960) the paths of the low-energy cosmic ray particles approaching the earth 
iB the solar merecions In this paper, we shall calculate the motion of these particles 
in more details and then examine their orbits in the earth’s magnetic field and the 
features of their precipitation on the earth. 

At first, as shown in Fig. 1, let the latitude of observin 
and the asymptotic velocity vectors of the cosmic ra 
upon the earth be (Ax, ¢.), where 4. and Poo 


g points be 4 (degrees), 
y particles impinging vertically 
are latitude and longitude components 
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Particle Orbit 


Fig. 1. Angle used to describe orbits of low energy 
cosmic rays coming from the sun and arriving vertically 
upon the earth, (from Firor (1954)) 

(a) Orbits in a meridian plane that moves with particle, 
showing the latitude of incidence, 4, the latitude of 
the asymptotic velocity vector, 20; 

(b) Projection of orbit on the equatorial plane. The 
angle gu is the longitude of impact relative to the source. 


of these velocity vectors, respectively. As valuable works on these velocity vectors 
of the cosmic ray particles within Bev-rigidity range were already performed by 
many authors (Malmfors, 1945; Dwight, 1950; Schlititer, 1951; 1958; Firor, 1954; 
Brunberg, 1956; Jory, 1956; Liist, 1957; 1958; Liist and Simpson, 1957), we shall here © 


Magnetic Rigidity in Ge. 
3 


Latitude of Incidence 
Fig. 2. Relation between the geomagnetic 
latitude of impact on the earth of particles 
arriving from the sun with magnetic rigidi- 
ties of 30 Mev to 3 Bev, and the geo- 
magnetic latitude of source (Sakurai, 1960). 


Rigidity 


30 Ney 
42 Bev 


0° 60° 20° 140° A4#0° 300° 360° 
too 
Fig. 3. Relation between the geomagnetic longitude 
of impact on the earth of particles arriving from 
the sun with magnetic rigidities of 30 Mev, 1.2 Bev, 
and 10 Bev, and the latitude of the source 
(Sakurai, 1960). 
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calculate the motion of the cosmic ray particles with low energies of Mev range. 
We shall show, in Figs. 2,3,4, and 5, the relation of (Aw, @.) of the cosmic ray particles 
arriving at 4 to these particle rigidities. 

The results about the particle rigidities 1, 2 and 3 Bev in Fig. 2 are due to Firor 
(1954). We see from Fig. 3 that the relation of 4. to g. does not depend upon the 
particle rigidities. This fact shows, if compared with Fig. 2., that, in similar manners 
as the incidence of cosmic rays with Bev-energies, these low-energy particles also 
approach the earth with some anisotropies from infinite distance. These impact 
zone effects may not appear very clearly on the observation since the features of 
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Fig. 4. Particle rigidities arriving at geomagnetic latitude 2 and the latitude component of the 
asymptotic velocity vectors, 2c. 
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eee 
their incidence quite depend on the local time of the onset of the incidence because 
of the discrepancy between the geomagnetic dipole axis and geographic one in the 
high latitude region. If we use the geomagnetic coordinate, however, those effects 
mentioned above will exist in the incidence of the low-energy solar cosmic ray 
particles to the earth. As some conclusions expected from those results are described 
in previous paper (Sakurai, 1960), we shall not here discuss any more about this 
problem. 

In the next place, we shall examine about the precipitation of the low-energy 
solar cosmic ray particles on the earth. We shall study the precipitation lines of 
their vertical incidence using the geomagnetic coordinates, and then show the results 
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Magnetic Rigi ad ty 


I. 30 Mev. 
IT. /00 Mev. 
G. 600Hev. 


Fig. 6. Precipitation lines of low rigidity particles 
coming from the sun on the earth. 
N: geomagnetic north pole. 


in Fig. 6. The center of this coordinate is the geomagnetic north pole, and we take 
the solar direction upward. These forms are the same as those of Stormer (1955), who 
used them to explain the formation of the auroral zone. This figure shows in what 
manner the low-rigidity particles coming near the earth along the straight paths in 
the solar direction will impinge upon the earth. As evident from this figure, their 
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Fig. 7. Comparisons of the calculations with the observations (I); September 13, 1957. 


© shows solar direction. 
(a) Precipitation lines of low-rigidity particles (30 Mev, 100 Mev and 600 Mev) 


(b) The incidence of these particles before the onset of geomagnetic storm (7 hours before) 
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Fig. 8. Comparisons of the calculations with the observations (II); February 11, 1958. 
(a) The same as Fig. 7 (a). 
(b) 11 hours before the onset of geomagnetic storm. 


Fig. 9. Comparisons of the calculations with the observations (III); February 11, 1958. 
(a) The same as Fig. 7 (a). 
(b) 6 hours before the onset of geomagnetic storm. 


precipitation lines are circular around the geomagnetic north pole. The greater the 
particle rigidities are, the larger become the radii of the precipitation lines. We 
shall show, in Figs. 7, 8 and 9, some examples of the comparisons of these theoretical 
results with the precipitation patterns obtained by the analyses of the observations, 
referring to the studies by Hakura, Takenoshita and Otsuki (1958) and Obayashi 
and Hakura (1960). As evident from the above comparisons, the incidence of the low- 
energy solar cosmic ray particles to the earth before the onset of geomagnetic Sc 
storms agrees very well with the theoretical results. Accordingly, we may see that 
the interplanetary space long before the onset of geomagnetic Sc storms is magnetically 
clean so as to be non-effective on the propagation of these particles. 

There are, however, some instances which are against what has just been said. 
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Since there are some observational cases where these low-energy cosmic ray particles 
impinge isotropically upon the earth even at the initial stage of their incidence (Fig. 
10), we cannot necessarily say that the interplanetary space are always magnetically 


ii 


Fig. 10. The incidence of low-energy solar cosmic ray particles on Oct. 21, 1957. This figure 
shows the isotropic incidence at the initial stage, and differs from Figs. 7, 8 and 9. 


clean. Consequently, such examples as shown in Figs. 7, 8 and 9 might be rather 
exceptional. Some recent studies on the low-energy solar cosmic ray events seem 
indeed to show that the incidence of these particles is isotropic from the initial stages 
(Reid and Leinbach, 1959; Leinbach, 1960), but, because the statistical analysis by 
Hakura (1960) on some low-energy solar cosmic ray events shows the clear existence 
of their anisotropic incidence before the onset of geomagnetic Sc storms, we may see 
that, statistically speaking, the low-energy solar cosmic ray particles will impinge 
upon the earth rather anisotropically at the initial stage of their incidence. 

Since the characteristics of the incidence of these low-energy cosmic ray particles 
to the earth are consistent with the observations on solar cosmic rays of Bev-energy 
range, the propagation mechanism of these low-energy solar cosmic ray particles may 
be also similar to that of the solar cosmic ray unusual increases so far observed 6 
times (Cocconi, Gold, Greisen, Hayakawa and Morrison, 1958; Obayashi and Hakura, 
1960; Sakurai, 1960). Hence, we may infer that the impact zone effects of the incidence 
of the low energy solar cosmic ray particles do not appear sharply because of the strong 
interaction in their passage with the interplanetary magnetic fields. 

We shall show in Fig. 11 a few examples of the precipitation lines of solar cosmic 
rays of Bev-energy range with respect to 1, 2 and 3 Bev particle rigidities at a fixed 
time. Fig. 1la shows the case of the solar cosmic ray event at the initial stage on 
February 23, 1956, while Fig. 11b shows the case on November 19, 1949. The hatched 
areas in these figures show the impact zones in both events mentioned above. 

We have examined, in this section, the characteristics of the incidence of the 
low-energy solar cosmic ray particles to the earth before the onset of geomagnetic 
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Fig. 11. The examples of the precipitation lines of Bev-rigidity solar cosmic 
rays at the initial stage. 
(a) February 23, 1956 (b) November 19, 1949 


Sc storms. And then, we could see that these particles impinge upon the earth in 
accordance with the expectation from the calculated results on their orbits (Figs. 7, 
8 and 9). On the other hand, there are some cases in which their orbits are much 
disturbed by the influence of the interplanetary magnetic fields (Fig. 10). We can, 
therefore, conclude that the propagation mechanism of these low-energy solar cosmic 
ray particles may always be controlled by the physical states of the interplanetary 
space. In the next section, we shall consider the behavior of these low-energy solar 
cosmic ray particles after the onset of geomagnetic Sc storms. 


3. Geomagnetic Storm Effects upon the Incidence of Low-Energy 
Solar Cosmic Ray Particles 


According to the studies of Hakura, Takenoshita and Otsuki (1958) and Obayashi 
and Hakura (1960), the incident region of the low-energy solar cosmic ray particles, 
after the onset of geomagnetic Sc storms, extends southwards beyond the polar cap 
region and shows spiral structures. Since the southward shifting of this incident region 
is very remarkable, it seems that there is a mechanism, in the outer space around the 
earth, quite different from that of the incidence of these particles before the onset 
of geomagnetic Sc storms discussed in section 2. 

We shall consider in this section the geomagnetic storm effects on such incidence 
of the low-energy solar cosmic ray particles. It is clear that the earth’s dipole 
magnetic field extends fairly far (say, 13 earth radii) if the geomagnetic storm 
producing solar plasma does not surround the earth (Sonett, Judge and Kelso, 1959; 
Sonett, Judge, Sims and Kelso, 1960). Such an extension of the earth’s dipole magnetic 
field was pointed out theoretically by Dungey (1955), and furthermore, the interplanetary 
plasma appears to confine the geomagnetic fields within the domain around the earth. 
It is, at present, inferred that, if the geomagnetic storm producing plasma surrounds 


Motion of Low-Energy Solar Cosmic Ray Particles in the Earth’s Magnetic Field 67 


the earth, the earth's magnetic fields will be confined within the domain of several 
earth radii from the earth (Chapman and Ferraro, 1932; Parker, 1958c; Dessler and 
Parker, 1959). 

In consequence, we can easily see that the distribution of the earth’s magnetic 
field in the earth’s outer space during the geomagnetic Sc storms may be quite 
different from that on the magnetically quiet conditions. We may, therefore, suppose 
that, if the low-energy solar cosmic ray particles interact with the earth’s magnetic 
field in its disturbed state, the conclusions in section 2 will be considerably altered. 
Hence, if these cosmic ray particles enter the severely disturbed earth’s magnetic 
field, their paths will be largely deflected, and their directions will become at random. 
We must, therefore, assume some mechanisms in order to explain the characteristics 
of their incidence during geomagnetic storms. Although such trials have been done 
in details by Ray (1956), Rothwell (1959), Obayashi and Hakura (1960) and others, 
this problem has not been solved satisfactorily. In order to make clear the mechanism 
of their incidence during the geomagnetic storms, we must explain both of the absence 
of their incidence in the polar-cap region and the cause of the spiral pattern, which 
is a remarkable feature in the incidence during geomagnetic storms (Hakura, 
Takenoshita and Otsuki, 1958; Obayashi and Hakura, 1960). For the complete 
explanation of this feature, it is necessary to know the states of the earth’s outer 
atmosphere in detail. 

Recently, Dungey (1958) considered theoretically the distribution of the magnetic 
fields of stars within their outer atmosphere and showed that these magnetic fields 
were confined within the restricted domain around them by the interaction with 
interstellar plasmas. As shown by Ferraro (1937), Liist and Schltiter (1955) and Cowling 
(1957), this result is also proved from the theoretical studies on the stellar magnetic 
fields. In view of these studies, we may infer that the forms in the outer space of 
the magnetic lines of force intersecting with the polar latitude region of the earth are 
quite different from those with low latitude regions (Hines, 1960). We may, therefore, 
suppose that very complicated phenomena will occur in the polar latitude region on 
geomagnetically disturbed conditions. And then, a configuration of the earth’s magnetic 
field, which forbids the incidence during geomagnetic storms of the low-energy solar 
cosmic ray particles to the polar latitude region, will be formed in the earth’s outer 
space. We do not yet know, at present, the theory by which the motion of the low- 
energy solar cosmic ray particles in the earth’s magnetic field during geomagnetic 
storms can be explained fairly. A qualitative consideration on this problem is presented 
by Winckler and Bhavser (1960), but, as this cannot explain it in details, we must 
continue to study the mechanism of their incidence during geomagnetic storms. 

In concluding, we must hereafter proceed to examine this unsolved problem since 
we can explain, as shown in section 2, the characteristics of their incidence before 
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4, Concluding Remarks 


The low-energy cosmic ray particles produced at the sun will propagate into the 
interplanetary space and then, interacting with the earth’s magnetic field, they will 
impinge upon the polar latitude region of the earth. We showed in section 2. that 
these cosmic ray particles impinge upon the earth at the initial stage of the incidence 
before the onset of geomagnetic Sc storm in accordance with the expectation from 
the theoretical results on their orbits (Figs. 7, 8 and 9). Hence, we can see that, 
when these particles propagate into the interplanetary space, they will approach the 
earth without the directions of their paths being much disturbed. Since we also know 
some cases where they have already impinged isotropically upon the earth even at 
the initial stage (Fig. 10), we cannot conclude that the interplanetary space is always 
transparent for their propagation. The propagation mechanism shown in a previous 
paper (Sakurai, 1960) will be able to explain the above observations. 

In consequence, the investigations on the low-energy solar cosmic ray events 
will be powerful means to make clear the properties and variabilities of the physical 
states of the interplanetary space. In particular, as shown in section 3, we do not 
stand on the satisfactory ground for the explanation of the incidence of these cosmic 
ray particles during geomagnetic storms, so we must continue to investigate this 
problem in future. The solutions of the mechanism of geomagnetic Sc storms and 
states of the earth’s outer atmosphere during these storms will settle this problem 
completely. 


Acknowledgement 


The author wishes to express his hearty thanks to Prof. Y. Tamura and Prof, M. 
Hasegawa for their kind advice and encouragement, and to Prof. K. Maeda and Dr. 
H. Maeda for their helpful discussions and criticisms on this work. He is also due to 


Mr. H. Leinbach, Dr. T. Obayashi and Dr. Y. Hakura for their valuable discussions 
and suggestions. 


References 


Anderson K.A. (1958) Phys. Rev. Lett. Li a35y 


Anderson K.A., Arnoldy R., Hoffman R., Peterson L. and Winckler J.R. (1959) J. Geophys. Res. 
64, 1133. 


Bailey D.K. (1957) J. Geophys. Res. 62, 431. 

Bailey D.K. (1959) Proc. IRE. 49, 255. 

Brunberg E. (1956) Tellus 8, 215. 

Chapman S. and Ferraro V.C.A. (1932) Terre. Mag. 37, 147. 

Chapman S. and Ferraro V.C.A. (1932) Terre. Mag. 37, 421. 

Cocconi G., Gold T., Greisen K., Hayakawa S. and Morrison P. (1958) Nuovo Cim. Supp. 8, 161. 
Cowling T.G. (1957) Magnetohydrodynamics ; Interscience Publ. 

Dessler D.J. and Parker E.N. (1959) J. Geophys. Res. 64, 2239. 

Dungey D.W. (1955) Rep. Phys. Soc. Conf. on Phys. of Ionosphere at Cavendish Lab. Cambridge, 


Motion of Low-Energy Solar Cosmic Ray Particles in the Earth’s Magnetic Field 69 


Sep. 1954. 

Dungey D.W. (1958) Cosmic Electrodynamics ; Cambridge Univ. Press. 

Dwight K. (1950) Phys. Rev. 78, 40. 

Ferraro V.C.A. (1937) M.N. 97, 458. 

Firor J. (1954) Phys. Rev. 94, 1017. 

Freier P.S., Ney E.P. and Winckler J.R. (1959) J. Geophys. Res. 64, 685. 

Gold T. (1959) J. Geophys. Res. 64, 1665. 

Gold T. (1960) Ap. J. Supp. 4, 406. 

Hakura Y., Takenoshita Y. and Otsuki T. (1958) Rep. Ionos. Res. Japan 12, 459. 

Hakura Y. and Goh T. (1959) J. Radio Res. Lab. 6, 635. 

Hakura Y. (1960) Lecture at the Society of Terre. hai and Electr. Japan, May 1960, Tokyo. 

Hines C.O. (1960) J. Geophys. Res. 65, 141. 

Jory F.S. (1956) Phys. Rev. 103, 1068. 

Leinbach H. (1960) Private Communication. 

Leinbach H. and Reid G.C. (1959) Phys. Rev. Lett. 2, 61. 

Liist R. (1957) Phys. Rev. 105, 1827. 

Liist R. (1958) Nuovo Cim. Supp. 8, 176. | 

Liist R. and Schliiter A. (1955) Z.f. Astrophys. 38, 190. 

Liist R. and Simpson J.A. (1957) Phys. Rev. 108, 1563. 

Malmfors K.G. (1945) Ark. for Mat. Astr. Fys. 32A, 1. 

Meyer P., Parker E.N. and Simpson J.A. (1956) Phys. Rev. 104, 768. 

Ney E.P., Winckler J.R. and Freier P.S. (1959) Phys. Rev. Lett. 3, 183. 

Obayashi T. and Hakura Y. (1960) J. Radio Res. Lab. 7, 27. 

Parker E.N. (1958a) Ap. J. 128, 664. 

Parker E.N. (1958b) The Plasma in a Magnetic Field, p. 77, ed. by Landshoff R.K.M.; Stanford 
Univ. Press. 

Parker E.N. (1958c) J. Phys. Fluids 1, 171. ‘ 

Ray E.C. (1956) Phys. Rev. 101, 1142. 

Reid G.C. and Leinbach H. (1959) J. Geophys. Res. 64, 1801. 

Rothwell P. (1959) J. Geophys. Res. 64, 2026. 

Sakurai K. (1960) J. Geomag. Geoel. 11, 152. 

Schliiter A. (1951) Z.f. Naturforsch. 6a, 613. 

Schliiter A. (1958) Nuovo Cim. Supp. 8, 349. 

Simpson J.A. (1960) Ap. J. Supp. 4, 378. 

Sonett C.P., Judge D.L. and Kelso J.M. (1959) J. Geophys. Res. 64, 941. 

Sonett C.P., Judge D.L., Sims A.R. and Kelso J.M. (1960) J. Geophys. Res. 65, 55. 

Stormer C. (1955) Polar Aurora; Oxford Univ. Press. 

Thompson A.R. and Maxwell A. (1960) Planet. Space Sci. 2, 104. 

Winckler J.R. (1960) J. Geophys. Res. 65, 1331. 

Winckler J.R. and Bhavser P.D. (1960) J. Geophys. Res. 65, 2637. 


On the Relativistic Electrons in the Solar Atmosphere 


By Kunitomo SAKURAI 


Geophysical Institute, Kyoto University 
(Read October 31, 1960; Received December 10, 1960) 


Abstract 


We consider the various processes, in solar atmosphere, of the relativistic 
electrons ejected in association with flares and examine how these electrons will 
lose their energies. 

In view of the above examination, we can show that these electrons will be 
trapped in the sunspot magnetic fields (~100 gauss) and then lose the majority 
of their energies through the synchrotron radiation. So, the radio waves emitted 
by these synchrotron radiation processes will be observed as the type IV outbursts. 


1. Introduction 


The behavior of the electron component associated with the acceleration of solar 
cosmic rays does not yet seem to have been made clear in details. 

Recently, however, some theoretical studies on the emission mechanism of the 
solar radio type IV outbursts have supplied some new informations of the motion of 
high energy electrons in solar atmosphere (Boischot, 1957, 1959; Boischot and Denisse, 
1957). At present, low-energy solar cosmic rays with energies ranging from some 10 
Mev to several 100 Mev are being investigated, and the relation of these cosmic 
rays to the solar and terrestrial disturbances, or the propagation mechanism of these 
cosmic rays into interplanetary space and so on, have been made clear both experi- 
mentally and theoretically (Reid and Leinbach, 1959; Obayashi and Hakura, 1960; 
Winckler, 1960). 

In this paper, we shall at first examine various processes, in solar atmosphere, 
of the electrons ejected in association with flares. And then, we shall consider the 
behavior of the relativistic electrons in solar atmosphere. Next, we shall present a 
possible explanation of the solar radio type IV outbursts. 


2. Relativistic Electrons in Solar Atmosphere 


The production of low-energy cosmic rays from the sun associated with flares 
was recently made clear by many authors (Bailey, 1957, 1959; Reid and Leinbach, 
1959 ; Simpson, 1960; Winckler, 1960). Furthermore, the production of low-energy solar 
cosmic rays is closely associated with the emission of the solar radio type IV out- 
bursts (Boischot and Denisse, 1957; Hakura and Goh, 1959). 


These facts suggest that the existence of the protons with energies ranging from 
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a few Mev to several 100 Mev in solar atmosphere will show that of the electrons 


with energies of the same orders. It will, therefore, be necessary to examine how 


the electron component will behave during the acceleration of solar cosmic rays. 
Both of protons and electrons ejected with solar flares will be accelerated or 
decelerated through their interaction with plasma and electromagnetic fields in solar 
atmosphere. Since the rates of energy loss of these protons through the interaction 
with solar plasma are negligibly small as shown by Hayakawa and Kitao (1956) and 
Parker (1957), it is clear that the acceleration of the proton component in solar 


atmosphere is possible enough. Then, for an acceleration mechanism, the Fermi 
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Fig. 1. Various processes of energy loss of the relativistic electron. 
The energy loss rates are shown with respect to electron kinetic energy. 
: synchrotron radiation loss. 

—--—-: ionization loss. 

------ : radiation loss. 
n, H, show the electron number density per cm? and magnetic field 
intensity (gauss) in a medium. 
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mechanism is proved to be most effective (Fermi, 1949; Parker, 1957, 1958a). Hence, 
the proton component will be well accelerated to cosmic ray energies. Furthermore, 
the studies on the suprathermal particles have recently made a progress, and their 
results are applied to explain the acceleration of low-energy solar cosmic rays and 
the observed characteristics of these particles (Parker, 1958b; Parker and Tidman, 
1958; Obayashi and Hakura, 1960). On the contrary, it is clear from Fig. 1 that the 
electron component will be decelerated rather than accelerated because its ionization 
energy loss rate is very large in solar atmosphere. Our following discussion, therefore, 
will be limited on the deceleration of these electrons. 

The electrons ejected into solar atmosphere in association with a flare will undergo 
various interactions with magnetic clouds there. Because the rate of electron accelera- 
tion is negligible, we shall examine various deceleration processes such as ionization 
energy loss, radiation effect and synchrotron radiation, and estimate these energy loss 
rates using electron number densities and magnetic field intensities as parameters. 
In order to estimate these rates, we here use the formulae derived by Heitler (1954), 
Ginzburg (1958) and Hayakawa, Ito and Terashima (1958). These estimations are shown 
in Fig. 1. From this figure we can see that, if the electron energies are in Mev-range, 
the energy loss rate by synchrotron radiation is most important in the strong magnetic 
field in outer atmosphere. In other words, as the electron densities in solar outer 
atmosphere (~a few solar radii) are very small, the energy loss rates by ionization 
and radiation become negligibly small and, in the strong magnetic fields (100 gauss), 
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Fig. 2. The distribution of the electron number 
density in the solar outer atmosphere. 
R@ shows the radius of the sun (=7X 10cm). 
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the synchrotron radiation loss process is more important than the above two processes. 
In the solar inner atmosphere, however, the ionization process is very effective as 
seen from Figs. 1 and 2. The data shown in Fig. 2 are due to Allen (1955). The 
radiation process is negligible in the outer atmosphere. In case of the magnetic field 
intensities are <10 gauss, the ionization process is most important and, if the intensity 
is greater than 100 gauss, the synchrotron radiation process becomes most effective. 
We note here that the electron kinetic energy range in the above discussion is from 
10 Kev to 1000 Kev (=1 Mev). 

From the previous discussion, we may infer that, after ejected into solar atmos- 
phere, the proton component will be accelerated to cosmic ray energies and propagate 
into interplanetary space and that the electron component will be decelerated mainly 
through the ionization process in the weak magnetic field regions or through the 
synchrotron radiation process in such strong magnetic field regions as sunspots (100 
gauss). 

| Since it has been guessed that the solar radio type IV outbursts are due to the 
synchrotron radiation of relativistic electrons (Boischot, 1957, 1959), we can see that 
the type IV outbursts are made by the electromagnetic radiation emitted through the 
interaction of the relativistic electrons with the sunspot magnetic fields because of the 
close association of the solar flares with the sunspot magnetic field configuration 
(McLean, 1959; Takakura, 1959). Pretty high correlation of the solar radio type IV 
outbursts with the low-energy solar cosmic ray events may show that both of protons 
and electrons will undergo the processes mentioned above (Hakura and Goh, 1959; 
Obayashi and Hakura, 1960; Sakurai, 1960; Thompson and Maxwell, 1960). 


3. A Possible Explanation of Type IV Outbursts 


As shown in the last section, the relativistic electrons ejected into solar atmosphere, 
roeiied with flares, will lose the majority of their energies within the domain of 
several solar radii from the sun while emitting the electromagnetic radiation. Thus, 
these relativistic electrons will disappear in solar atmosphere, so they will not propagate 
into interplanetary space. The large dispersion of the emitting sources of type IV 
outbursts, ranging from 1.2 to about 5 solar radii, leads us to conclude that the emit- 
ting sources will move actively in solar atmosphere. And then, as the number of 
relativistic electrons ejected from the sun will be extremely large in the low energy 
end, considerably strong magnetic field of the order of 100 gauss must exist in the 
solar outer atmosphere in order to explain the observations of the type IV outbursts. 
A possible explanation, therefore, is as follows: After the high energy protons and 
electrons are ejected into solar atmosphere in association with flares, they will transport 
upward the surrounding sunspot magnetic fields. Thus, we may suppose that, even 
within the solar outer atmosphere, the strong magnetic field of several 100 gauss 
exists at the time of these disturbances. These electrons will consequently lose their 
energies through the synchrotron radiation which is caused by their interaction with 
the sunspot magnetic fields stretched out by the solar plasma. So, we may observe 
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Fig. 3. The ratio of the proton magnetic rigidity to the electron’s 
shown with respect to the kinetic energy (see Appendix 1). 


these processes as solar radio type IV outburst. 

In Fig. 3, we shall show the ratio of the proton magnetic rigidity to that of 
electron in regard to the same kinetic energy (Appendix 1.). Since this ratio is very 
large as seen from this figure, it seems that the trapping of these electrons in magnetic 
fields is much easier than that of the proton. Hence, we may interpret that, while 
the proton component will escape more easily from the magnetic fields extended into 
the solar outer atmosphere, the electron component will be unable to do so and they 
will lose there the majority of their energies. 

The duration, so-called half-life, of the emissions of the type IV outbursts is 
ranging from about 20 minutes to 1.5 hours. Taking into account these values, we 
shall show in Fig. 4 the relation of the electron initial kinetic energies to the magnetic 
field intensities H, (gauss unit), where H, is a component of magnetic field (H) 
perpendicular to the electron velocity vector. We can also see from this figure that 
the relativistic electrons will lose their energies mainly through the synchrotron radia- 
tion process. 

In concluding, we can see that the cause of the type IV outbursts is due to the 
interactions of the relativistic electrons of energies of some ten to several hundred 


Kevs, associated with solar flares, with the sunspot magnetic fields extended into the 
outer atmosphere. 


4. Conclusion 


After ejected into solar atmosphere in association with solar flares, component of 
both protons and electrons will diffuse into the outer space, undergoing various processes. 
The proton components will propagate into interplanetary space after accelerated to 
cosmic ray energies. On the other hand, since the deceleration processes are very 
important for the electron components, they will lose their energies through the 
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Fig. 4. Relation of the magnetic field intensity (H,) with 
the electron initial kinetic energy, where the half-life 
of synchrotron radiation (714) are used as parameters. 


ionization and the synchrotron radiation. 

We can interpret the observations of the type IV outbursts on the basis of the 
assumption that there are strong magnetic fields of the order of 100 gauss in the 
solar outer atmosphere, which are extended by solar plasmas; in other words, these 
outbursts will be produced by the interaction of electrons with the sunspot magnetic 
fields. 


Acknowledgement 


The author wishes to express his sincere thanks to Prof. M. Hasegawa and Prof. 
Y. Tamura for their kind advice and encouragement. His thanks are also due to 
Prof. K. Maeda and Dr. H. Maeda for their valuable discussions and criticisms through- 
out the study. 


Appendix 1. 


Ratios of the Proton-to the Electron Magnetic Rigidity 


Kinetic energy, momentum and magnetic rigidity of a charged particle (charge Ze; 


rest mass m) are expressed as follows: 


1 
kinetic energy Eug=me| ae —1| ; (1) 


76 K. SAKURAI 


momentum P= itis ; (2) 
. segs ti Ein Exint+2mc?) hae 
and magnetic rigidity Fe =|. Evin (Ze)? — } , (3) 


where c is light velocity and $= > v being a particle velocity. 


We shall use here a suffix “pf” and “e” for proton and electron respectively. 
The ratio of the proton rigidity to that of electron can be shown in the following 
form, if kinetic energies of both particles are quivalent (Z7=1). 

Proton magnetic rigidity =| Fewt inne)", (4) 
Electron magnetic rigidity Ein t+ 2mec* 


where Ejinsp=Egine=Lxin: 

The results deduced from equation (4) are shown in Fig. 3. We can see that, 
as the electron rigidity is much smaller than the proton’s even if the electron kinetic 
energy is the same as the proton’s, the electrons will be easily trapped in the magnetic 
fields. 
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Abstract 


At present, there are some theories being able to explain the frequency 
characteristic of dawn chorus. In these theories, the natural thermal noise 
excited in a manner similar to the operation of a travelling-wave tube, the 
Doppler shifted proton cyclotron radiation and the Cerenkov radiation are considered 
as the origin of V.L.F. noise in the earth’s exosphere respectively. This paper 
indicates that the Cerenkov radiation will be the most effective with respect to 
the emission energy in the above mentioned mechanisms in the earth’s exosphere 
and that Cerenkov radio emissions generated by high speed protons (1.4 x 10 km/sec) 
precipitating into the ionized exosphere along the earth’s magnetic field may be 
excited simultaneously by the operation of the mechanism similar to a travelling- 
wave tube due to charged particles different from those causing the Cerenkov 
radiation relatively near the earth less than about 4 earth radii from the center 
of the earth and then propagate toward the earth in the whistler mode. 


1 Introduction 


Storey (1953) has pointed out that the audio- frequency phenomenon .called “ dawn 
chorus” is the sound of a distant origin and that its occurrence correlates strongly 
with geomagnetic activity. 

Such a same result has been also reported by Allcock (1957). 

Then, latitude effect on the diurnal maximum of dawn chorus activity has been 
found by Allcock (1957) and Pope (1957, 1960) as represented in Fig. 1. Allcock has 
made the hypothesis that dawn chorus signals are generated in the vicinity of the 
geomagnetic equatorial plane. 

He has replotted the data of Fig. 1 as a polar graph of local geomagnetic time 
against the height above the equator of the geomagnetic line of force passing through 
the station and has shown that the curve fitting the data is characteristic of the locus 
of a positively charged incoming particle in the equatorial plane, and being deflected 
by the earth’s magnetic field. Thus, from above mentioned facts we can certainly 
deduce that dawn chorus may be generated directly or indirectly by streams and 
bunches of high speed charged particles precipitating into the ionized exosphere in 
the presence of the earth’s magnetic field. Next we review briefly some theoretical 
works concerning the emission mechanism of V.L.F. noise in the exosphere. 

Gallet and Helliwell (1956, 1959) have proposed that V.L.F. noise emissions are 
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due to the selective amplification of low-level electromagnetic waves, such as caused 
by the natural thermal noise in the exosphere. This amplification mechanism is similar 
to the operation of a travelling-wave tube and the condition for the operation of this 
mechanism is completely satisfied relatively near the earth, less than about 4 earth 
radii. 

Recently Barrington (1959) has suggested that the condition for coupling between 
the whistler mode and the space charge modes of the beam resulting in an amplification 
of the whistler mode seems likely for certain region of the exosphere, especially above 
2000 km over the ground. 

Thus it is appropriate that we seek the region, where the interaction between 
the whistler mode and the beam modes takes place, relatively near the earth above the 
ionosphere. MacArthur (1959) has suggested that dawn chorus may be due to Doppler- 
shifted radiation from protons which are rotating about the geomagnetic lines of 
force with the cyclotron frequency, but he has argued it concerning the frequency 
characteristic only. Veksler (1957) has suggested that the Cerenkov radio emissions are 
generated by charged particles moving along magnetic lines of force in plasma. It 
has been recently indicated by Ellis (1959) that V.L.F. phenomena are due to the 
Cerenkov radiation by high speed charged particles moving along the geomagnetic 
lines of force. 

Here we discuss the effectiveness of above mentioned emission mechanism in the 
earth’s exosphere from the point of view of emission energy and compare the emission 
energy of a possible more effective mechanism with observed values of V.L.F. noise 
at ground. 


LATITUDE VARIATION IN TIME OF DIURNAL 
MAXIMUM OF DAWN CHORUS ACTIVITY 


ACTIVITY (LOCAL GEOMAGNETIC TIME) 


TIME OF MAXIMUM DAWN CHORUS 


50 60° 
GEOMAGNETIC LATITUDE 


70° 


Fi 


— 


g. 1. Latitude Variation in Time of Diurnal Maximum 
of Dawn Chorus Activity (after Allcock) 


2 The Origin of Dawn Chorus Deduced from Observed Facts 


Allcock (1957) has found a marked variation of the time of maximum activity 
of dawn chorus with geomagnetic latitude from many data as shown in Fig. 1. Now 
we make the hypothesis that dawn chorus signals are generated relatively near the 
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Fig. 2. 


earth above the ionosphere by charged particles penetrating into the earth’s exosphere. 
It is then appropriate to replot the data of Fig. 1 as a polar graph of local geomagnetic 
time of the diurnal maximum of dawn chorus activity against the geomagnetic latitude 
of each station. This is done in Fig. 2 with open points and a dotted curve. We 
can see that a dotted curve in Fig. 2 represents a spiral expanding in a clockwise. 

On the other hand, Meek (1955) has used values of hourly averages for the year 
of “disturbed days minus quiet days” which Vestine (1947) has tabulated and charted 
from the 1932-33 Polar-Year data, and has made a ‘polar plot of the local geomagnetic 
time of the diurnal maximum decrease in geomagnetic horizontal component for all 
available magnetic stations north of 40°N. This is represented by solid points and a 
shaded band that is also a spiral expanding in a clockwise in Fig. 2. These spirals 
expanding in a clockwise are characteristic of the Stérmer’s precipitation spiral of 
positively charged particles and we can see that the spiral concerning dawn chorus 
shifts toward lower latitudes than the spiral concerning geomagnetic horizontal com- 
ponent. 

It is well known that geomagnetic disturbances at high latitudes are caused by 
charged particles impinging upon the ionosphere and the auroral zone gradually shifts 
southward as the corresponding storm grows intense (Nagata, 1949). From above 
facts, we can deduce that dawn chorus signals are generated by considerable high 


energy positively charged particles in the solar corpuscular stream. 


3 A Possible More Effective Emission Mechanism of V.L.F. 
Noise in the Earth’s Exosphere 


Now we discuss which the proton cyclotron radiation or the Cerenkov radiation 
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due to protons is more effective with respect to the emission energy in the exosphere. 
We consider the longitudinal extraordinary propagation of V.L.F. noise in both cases 
in the exosphere. 

The emission energy of the Cerenkov radiation and the proton cyclotron radiation 
at the frequency of 5kc/s are estimated as 1.7X 10-24 watts. m=? and 7.2x10-*8 watts. 
m-® respectively at 2 earth radii from the center of the earth on the geomagnetic line 
of force passing through the earth’s surface at the geomagnetic latitude of 60°. About 
this point in the exosphere, the electron density is 5X10%/cm*, which is deduced from 
the extrapolation of whistler observation (Ondoh and Hashizume, 1960), the velocity of 
protons is 10®cm/sec, and the density of corresponding protons is 10-!/cm® which is 
of the order of auroral protons. 

Thus the emission energy of the Cerenkov radiation is much greater than that 
of the proton cyclotron radiation under the same conditions in the exosphere. Hereafter 
we shall consider only the Cerenkov radiation as the origin of V.L.F. noise in the 
exosphere. 

Then we review briefly the amplification mechanism similar to the operation of a 
travelling-wave tube given by Helliwell and Gallet (1956). By analogy with the small 
signal theory of travelling-wave tubes the necessary condition for coupling between 
the electromagnetic wave and the stream of charged particles is assumed to be V=vp 
where V is the velocity of the stream of ionized particles and vp phase velocity of 
the electromagnetic wave. 

For small propagation angle the quasi-longitudinal approximation holds good. 


In the whistler mode the local refractive index is n=[1+ 37/2" and vp=— 


where f is the wave frequency, f) the plasma frequency and /;, the electron gyro- 
frequency. 
Combining vp and », we obtain ca solutions of f which are given by 


bshi-(efe EP] 


For interaction to occur, f must be teal and hence the condition for interaction 


to occur is (240.2 \'< 1. The physical conditions of the exosphere for (2fet)"< 1 will 
be satisfied relatively near the earth less than about 4 earth radii from the center of 
the earth. 

eens note especially the condition for interaction to occur and the frequency 
characteristic. As Ellis (1959) has suggested, the Cerenkov radio emissions due to 
charged particles moving along the geomagnetic line of force give the same type 
frequency characteristic as that of the mechanism given by Gallet and Helliwell. The 


condition for the Cerenkov radio emissions to occur is given by Vcos@=& where V 
n 


is the velocity of charged particles moving along the geomagnetic line of force and @ 
the auzle between the wave normal direction and the geomagnetic line of force, c the 
velocity of light, and n the local refractive index for the extraordinary mode 
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+3 eG. = Fo 1/2 ; 
Combining Vcosd=—) and a | , We obtain the frequency 


of the Cerenkov radio emission of the extraordinary mode due to high speed 
charged particles moving along the Bi ecas line of force as 


f= yf xcosd] 14{1 —(24 oto 0.7 : Veal 
ce 


where /, is the plasma frequency, and Ff the electron gyrofrequency. 

see it 
och 1 —Vcos%a 
to unity. For the Cerenkov radio emissions of the extraordinary mode to occur, f must 
be real and therefore the condition for the generation of Cerenkov radio emissions is 


(242 ~\'< 1. 


H C 


Generally the value of the expression will be approximately equal 


As shown by Gallet, at distances less than about 4 earth radii from the center of 
4 2 - 
the earth, (242 =) will generally be very small compared to unity and so the 


Cerenkov radio emissions will occur. 
Under such conditions we adopt the lower frequency of solutions as f= fr e. ; i cosé 
HT 


(@3-0). In the anisotropic and dispersive medium, the angle between the energy flux 
direction of the electromagnetic wave and the geomagnetic line of force is generally 


smaller than that between the wave normal direction and the geomagnetic line of 


force in the whistler mode (Maeda and Kimura, 1956), and for large values of ( fo )" 


Shu 


that will be generally satisfied in the exosphere, energy flux directions are entirely 


limited to a cone of semi-angle 19°29’ about the field direction in the whistler mode 
(Storey 1953). ay 

Therefore the energy of the electromagnetic wave tends to travel in the direction 
of the particle stream as a result of the guiding action of the geomagnetic field. The 
condition of the generation of the Cerenkov radiation and the operation of amplifi- 
cation mechanism similar to a travelling-wave tube is given by V,cosd=~- and V,= 7" 
respectively in the medium of the local refractive index m, where V, and V2 are 
velocities of protons. 

When @ is small but not zero and m is greater than unity, we may consider that 
Cerenkov radio emissions of the whistler mode are amplified almost simultaneously 
by the operation of the mechanism similar to a travelling-wave tube due to charged 
particles different from those causing the Cerenkov radio emissions at different points 
in space and time. There is lateral stability of the beam when the Maxwell stress 
tensors of the geomagnetic field exceed the stress tensors of the self-fields and of the 
beam itself. Neglecting the self-field, this requires that H?/4n>Nm,V? where JN is 
the number density and m, the mass of an ion in the beam, and V is the velocity 
component parallel to the geomagnetic field H (Harrison 1960). 

In the part less than about 4 earth radii from the center of the earth on the 
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geomagnetic line of force passing through the earth’s surface at the geomagnetic 
latitude of 60°, the above condition is satisfied for auroral particles. 

Recently Neufeld (1959) has suggested that an electron beam passing through a 
dielectric medium may produce an instability that is associated with the growth of 
longitudinal waves having a velocity close to the velocity of beam and if inhomogenei- 
ties are present, these longitudinal waves may be converted into transverse waves and 
radiated into space at frequencies different from the Cerenkov frequencies, and also 
Barrington (1959) has suggested the possibility of the coupling between the whistler 
mode and the space charge modes of an ion stream in the exospheric region considered 
here. Now assuming that both the Cerenkov radio emissions and the amplification 
of these waves in a manner similar to the operation of a travelling-wave tube can 
occur in the exospheric region considered here, we calculate the energy of the Cerenkov 
radiation due to high speed protons at 2 earth radii from the center of the earth on 
the geomagnetic line of force passing through the earth’s surface at the geomagnetic 
latitude of 60° where is the most suitable part for above mentioned mechanisms to 
operate. 


The condition for the generation of the Cerenkov radio emissions (V><) is 


satisfied at above mentioned part where the electron density is 5x 10°/cm*, V=1.4x10* 
km/sec, and »=22.5 for the radio wave of the whistler mode of which frequency is 
5 kc/sec. We have for the total emission in a column of the length of ZL cm which 

4r*e? (1- 1 ) SYA 10# 
ce fPn® 10? 
where JN is the number density of protons, V the velocity of protons and pat (Shiff, 
1949, Marshall, 1956, and Ellis, 1957). 


The length of the column is limited within several hundreds kilometers by the 


consists of the high speed protons S= 


watts. m~*. (c/s)-! 


frequency bandwidth observed at a moment (Gallet, 1959). Using (er =I, Vs 


1.4X10°cm/sec, L=2x10’cm, N=10-!/cm*, f=5x10%c/sec, e=4.8X10- es.u., and 
c=3X10! cm/sec, we have S=1.4x 10-22 watts. m=*. (c/s)71. 

On the other hand, the incident power flux of dawn chorus recieved at the ground 
is of the order of 10-'® to 10-'* watts. m=?. (¢/s)"!. 

Considering that the Cerenkov radio emissions generated by incoming high speed 
protons are amplified in a manner similar to the operation of a travelling-wave tube 
and propagate in the whistler mode without a little loss on the way, we can not obtain 
theoretically the observed power flux of dawn chorus without a gain of amplification 
greater than 60 db. For the physical states existing in the outer ionosphere, in addition 
to condition of the matching of stream and wave velocities, Bell (1959) has obtained 


a gain of approximately 2 db per wavelength at frequency of 3kc/sec. This value 
of gain is very much smaller than we have required. 


Concluding Remarks 


As we have seen in latitude effect on diurnal maximum of dawn chorus activity, 
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it is clear that dawn chorus is generated in the earth’s exosphere by positively charged 
particles incoming from the outside of the earth. The most effective radio emissions 
at frequencies of V.L.F. noise in the exosphere may be the Cerenkov radiation due to 
high speed protons rather than the proton cyclotron radiation. In order to explain 
the V.L.F. phenomena, we shall have to consider the generation mechanism of V.L.F. 
noises and the amplification mechanism of these emissions at the same time in the 
exosphere. Ellis (1960) has recently reported that bursts of the V.L.F. noise occur even 
in geomagnetic undisturbed conditions. 

It seems likely that hydromagnetic mechanism to accelerate charged particles up 
to the order of the velocity of 10¢km/sec and unknown mechanism which frequently 
causes V.L.F. noises, exsist in the earth’s exosphere. 
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Abstract 


The isoporic charts of X, Y and Z of geomagnetic secular variation during 
the period from 1955 to 1960 are constructed based on data observed at 99 stations 
over the earth’s surface including the Antarctic region, and illustrated in Figs. 1-9. 
Equivalent current system of the variation on the earth’s core’s surface is shown 
in Figs. 11-13. 

From the results of the spherical harmonic analyses of these isoporic charts, 
it is concluded that (a) the earth’s magnetic dipole is decreasing in moment 
with the rate of —4.4x1022emu/yr, magnitude of the moment at 1960.0 being 
8.046 x 10? emu ; (b) the eccentric dipole is drifting westwards by 18’/yr, northwards 
by 12’/yr, and outwards by 3.4x10-*x the earth's radius/yr. It is shown that 
the other harmonics are drifting mostly westwards also with the rate of several 
decimals of a degree per year. 

On the other hand, very intense foci of isopors are found in the Antarctic 
region, amounting to about 200;7/yr in Z. The current change on the earth’s 
core’s surface equivalent to the foci amount to about 0.1 Amp/cm/yr. This value 
is strikingly large compared with intensity of the equivalent current of the earth’s 


magnetic dipole, which is about 0.5 Ampere/cm in this region. 


1. Introduction 


The general features of geomagnetic secular variation during the period from 
1950 to 1955 have already been reported on by one of the authors (T.N.) and T. 
Rikitake™. In continuation of the previous report, the general features of geomagnetic 
secular variation after 1955 will be summarized and analysed in the present report, 


together with some remarks on the results. 


At the request of one of the authors (T.N.), 87 standing magnetic observatories 
(excluding those in the Antarctic Continent and on some Sub-Antarctic islands) kindly 
sent him the annual mean values of the geomagnetic three components. In addition 
to the data from these established magnetic observatories, various information about 
geomagnetic secular variation in the Antarctic area could be obtained for the period 
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concerned. It seems that these Antarctic data are quite characteristic in this paper, 
because geomagnetic isoporic charts covering the whole parts of the earth’s surface 
could be compiled, based on these new data together with those in the other parts 
of the world. 

Thus, the total number of geomagnetic stations, whose data were used in the 
present work, amounts to 99. The positions of these magnetic stations are listed in 
Table I, where the periods of observation at each station refered to in the present 
work are also given, and are shown by full circles in Figs. 3, 8 and 9. 


2. Data 


At all standing magnetic observatories, including Macquarie Island, Orcadas del 
Sur and Mawson in the Antarctic region, annual mean values of geomagnetic three 
components were obtainable. For the purpose of the present work, where the variation 
during 1955-60 is dealt with, plots of the annual mean values against time at all 
observatories are classified into three categories, in which the annual rate of the 
variation has the following tendencies,-namely; (a) the annual rate has remained 
almost constant since 1950 or since the year before; (b) the magnitude of the rate 
changed more or less around 1955 and has kept nearly constant since then, but the 
tendency of the variation after 1955 still remains similar to that before 1955; and (c) 
the annual rate changed its sign after 1955; that is, the plotted curve has a maximum 
or minimum during the period concerned. Almost all data dealt with here belong to 
either (a) or (b) category, and only at six stations one of the three components can 
be classified in (c) category. In Table I, the data of category (b) are marked with a 
single asterisk (*), and those of category (c) with double asterisks (**), while those of 
category (a) with no mark. — 

As for data on categories (a) and (b), there may be no doubts in evaluating the 
mean annual rate of the variation for the period 1955-60, but there may be some 
ambiguity in estimating the mean rate for category (c) data. In the present work, 
however, the simple average of annual rates during the period concerned is adopted 
as the mean annual rate. 

On the other hand, most of the Antarctic stations were operated during IGY or a 
few years before and after IGY. By examining the monthly mean values, or in some 
cases the series of the results of absolute measurements, at these Antarctic stations, 
a fairly reliable value on geomagnetic secular variation could be estimated at 12 stations 
listed at the bottom of Table I. There are several other stations whose data are not 
sufficient for determining the secular variation, owing either to the too large effect of 
polar geomagnetic disturbances or to too short a period of observation. 

Since the IGY period was from July 1957 to December 1958, it seems that these 
data for the IGY period and its extension before and after can show approximately 
the annual rate of secular variation near the middle period of 1955-60, namely, 1957.5. 
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Table I. Annual Rate of Secular Variation in the Geomagnetic 
Field (X, Y, Z) during the Period 1955~60. 
Annual rate of 
CSAGI Geographic Secular Variation : 
sia Observatory Late Done. y Z Period of used data 
(y/year) 
° ° 1947~49 
Ao2l | Thule 76.5 290.9 1 ise 46 | Hostag 4 
A030 Resolute Bay fiingie <piasaal 29 iL 33 1954~55~58 
A031 Bjérndya 74.5 19.0 g* 20 39 1951~55~58 
A039 Point Barrow fie Se Usies' 3 —4 —-2 1950~55~59 
A047 Troms¢ 69.7 19.0 8 12 42 1945~55~58 
A049 Godhavn 69.2 306.5 19 7 41 1946~55~58 
A054 Abisco 68.4 18.8 4* 24 32 1946~50~56 
A092 College 64.9") “21222 133 —4 19* 1948~55~59 
A099 Baker Lake 64.3 3263-9 30 9 49* 1951~55~59 
Al21 Srednikan 62.4 152.3 8- —13 17 1946~55~58 
A123 Dombas 62.1 9.1 10 20 28 1955~58**** 
A124 Yakutsk 62.0 129.7 23 —21 13** | 1946~55~59 
A134 Nurmijarvi 60.5 24.7 | 16 33 1953~55~58 
A140 Lerwick 60.1 358.8 17 19 29 1946~55~59 
B002 Leningrad 60.0 3027 ae a 10 34 1948~55~59 
B009 Lovo 59.4 17.8 4 18 34 1949~55~58 
A149 Sitka 5720" 224270 17 — 7 — 3 1946~55~59 
BO19 Sverdlovsk 56.7 61.1 6-5. Sa DB 1946~55~59 
B114 Rude Skov 55.8 12.4 9 22 28 1946~55~59 
B028 Kazan 55.8 48.8 ae 5** 41 1946~55~58 
B035 Moscow 55.5 37.3 — 2 6 40 1946~55~59 
B038 Eskdalemuir 5b) o00.9 25 27 23 1946~55~59 
A154 Meanook 54.6 246.7 Sas ee 1950~54*** 
B058 Wingst 53.7 S51) 12 27 33 1945~55~59 
BO071 Witteveen 52.8 6.7 20 27 24 1945~55~59 
C362 Irkutsk 52.5 104.0 14 — 21 a 1946~55~59 
B089 Swider ya | bol begs! 1 21 48 1955~58 
B095 Niemegk Ly al too 10 25 32 1945~55~59 
B098 Valentia 51,95 73497 38 25 17 1945~55~59 
B14 Abinger 5Tj24 309.06 27 33 20 1946~55~59 1) 
B119 Hartland 5150). $55.5 32 vas 19 1957~59 
B132 Manhay 50.3 5.7 17 30 32 1951~55~59 
B136 Dourbes 50.1 4.6 21 24 28 1957~59 
B143 Pruhonice 50.0 14.6 11 27 37 1946~55~59 
B163 Fiirstenfeldbruck 48.2 ee 14 32 32 1946~55~58 
B172 Hurbanova 47.9 18.2 6 25 40 1949~55~ 
C016 Yugno-Sakhalinsk 46.9 142.7 1 .—10 —18* patina 
C018 Odessa 46.8 30.9 — 7* 14 34 1948~55~59 
B239 Castellaccio 44,4 8.9 19 38 25 1947~55~59 2) 
C034 Memambetsu 43.9 144.2 3 —9 —16**|} 1952~55~59 
B249 Agincourt 43.8 280.7 41 10 2 19 
‘ : . - 50~55~58 
oy Vladivostok 43.7" 181.9 TO a ae 1952~55~5713) 
ace Tbilisi 42.1 44.7 —]2*. — ]* 39 1946~55~59 
Tashkent 41.4 69.2 10 —1I* 35 1946~55~59 
Istanbul 41.1 29.1 O* 18 61 1949~55~59 
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Table I. (continued) 
BAC: a rate ut e 
I . te 
ins Observatory tee sae ae nae Period of used data 
(y/year) 
C093 Coimbra 40.2 351.6 43 43 — 8* | 1946~55~59 
C098 Toledo 39.9 355.9 35 46 —3 1947~55~59 
B314 Cheltenham 38.7 283.2 ee a 1946~50~56 
B318 Fredericksburg 38.2 282.6 43 -—-19 —34 1956~59 
San Migule 37.8 334.3 67 30 —39 1947~55~59 
C137 Almeria 36.9 357.5 42 56 —5 1955~58 
C143 San Fernando 36.5 353.8 44 52 1951~55~59 
C147 Kakioka 36.2 140.2 5 —7 —28* | 1946~55~59 
Ksara 33.8 35.9 12 11 63 1945~55~58 
C214 Simosato 33.6 135.9 6 —5 -5 1955~58 
C236 | Tucson 32.2 249.2 | -—10 -—11 —39 | 1946~55~59 
C251 Quetta _ 30.2 66.9 9 —16** 19 1954~55~59 
C256 | Helwan | 229 31.6 25 wm 48 { See 
C273 Tamanrasset 22.8 5.5 32 56 —11 1950~56 
C277 Honolulu -21.3 201.9 | -23 —4 —15* | 1946~55~59 
C287 Teoloyucan 19-8 260.8 == EL I5 10 1951 ~53*** 
Alibag 18.6 72.9 eh 8* | 1947~55~58 
C300 San Juan | 18.4 293.9 2* —A49 —103 1946~55~59 
Coil M’Bour 14.4 343.0 | 37 6l. —67 1952~55~59 
E553 Muntinlupa 14.4 121.0 | 136-217. =20 1951~55~59 
E556 Guam 13.5 144.8 3 6  =%6 1957~59 
E566 Kodaikanal 10> .. 775 GS 1950~55~59 
E575 Paramaribo 5.8 304.8 2] —48 +106 1958~59 
E583 | Bangui 46 6.6) en FPA eg 
E590 Tatuoca = 1.2) 311.5 fh =a7 +137 1958~59 
E625 Hollandia - 2.6 14.5 | —12 oe 1957~59 
Kuyper —60 106.7 | —6 4** _ 3 1954~55~59 
E644 Elisabeth ville +i 7) 2ES —35 105 =15 1945~50~55 
E646 Huancayo 12 Ne. 2847. —5A BG) — 3 1948~55~58 
E653 Apia —13.8 188.2 | —20* ll 6 | 1950~55~59 
C901 Tananarive -18.9 47.5 | —25 —18 11 1945~55~59 
Mauritius —201 57.6 | —20 -—23 -—20 | 1945~55~594) 
E702 La Quiaca —22.1 294.4 —9 —57 —9 1949~50~56 
E629 Vassouras ~22:4. 316.4 | —80. —38 . ~77 1946~55~58 
C925 Watheroo —30.3 115.9 0 6 —61 1951~55~589) 
C933 Gnangara —31.9 116.0 —4 -—-14 —36 1958~59 
C932 Pilar —31.7 296.1 264 482 8 1949~50~56 
C957 Hermanus —34.4 19.2 72, 26 65 1946~55~59 
B966 Toolangi —37.5 145.5 | —19 30 —5 | 1949~50~56 
B979 | Amberley | -43:2 172.7 | -28 24 19 | 1950~55~59 
A961 Macquarie Is. —54.5 159.0 —18 47 16 1952~55~59 
A962 Geman del Sur | —60.7 315.1 a6, 93 | 1933~55 
(Antarctic Region) \(1952~54~60 (for X & Y) 
Heard Is. 253.0 A | +37 —44 98 { 1952 ~54 (for Z) 
Ushuaia =§408 7251.7 | =—64. —29"" “120 1956~58 
Redeemed 63.4 303.0 | —57. —26 107 | 1956~58 
Merchior —64.3 297.0 | -49 -30 132 | 1956~58 
A973 Argentine Is. nS. 29BFi) 84, 48) 131 | 1956~58 
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Table I. (continued) 


Annual rate of 
G hic Secular Variation . 
Oe Cae Observatory eae 0% x y Z Period of used data 
(7/year) 

A977 Wilkes —66.3 110.5 —87 —1 21 IGY 
A978 Mirny —66.6 93.0 —94 —39 45 1956~58 
A980 Mawson —67.6 62.9 —71 —57 51 1955~58 
A984 Syowa —69.0 39.6 —79 —40 181 1957~59 
A989 Halley Bay —75.5 333.4 —13 —39 70 IGY® 
A991 Scott Base —77.9 166.8 —38 15 36 IGY 
A995 Little America —78.3 199.8 7 45 132 IGY 


Remarks: * Category (b) in § 2 of the text. 
** Category (c) in § 2 of the text. 
#k ~“Values before 1955 (for reference). 
#kke Average for quiet days only. 
(1) At Abinger until 1956, and at Hartland since 1957. 
Data obtained at Hartland since 1957 are reduced to those at Abinger. 
(2) Original data are H, D and J. 
(3) Data of 1958 and 1959 are ignored, because of their unreasonable disconti- 
nuity. 
(4) With fair amount of fluctuation on account of some trouble in magneto- 
graphs. 
(5) At Watheroo until 1958, and at Gnangara since 1958. 
(6) This station is built on ice shelf, which is believed to be moving with 
velocity of 400~600 m/year. Therefore the secular variation values are 
a little ambiguous. 


3. Isoporic Charts xX Y and Z for the Period 1955-60 


From the above-mentioned data, the mean annual rate of secular variation of the 
geomagnetic three components X, Y and Z were computed, the results being given 
in Table I and plotted in Figs. 1~9. In drawing isoporic lines for X, Y and Z, two 
basic assumptions were considered ; they are (A) that F =(X, y, z ) satisfies the con- 
ditions, curl F=0 and div F=0 near the earth’s surface, and consequently p2F'=0; 
and (B) that the origin of F is situated beneath the earth’s surface. In the practical 
procedure, an isoporic chart for Z was first constructed from the distribution of Z 
values with reference to the horizontal vectors H=(X, Y), based on the assumptions 
that div F=0 and (B), and then isopors for X and Y were drawn by considering 
0X/dy=0Y/dx. The result of the second process reflects more or less upon the Z 
contours drawn by the first process. Thus, by repeating successively the above-mentioned 
processes, isoporic charts for X¥, Y and Z Satisfying approximately the assumed con- 
ditions (A) and (B) were obtained. 

On examining Figs. 1-9 in comparison with the corresponding isoporic charts for 
1950-55, constructed by Nagata and Rikitake™, and also with those for 1940-45 by 
Vestine et al®, it is seen that the general tendency of geomagnetic secular variation 
over the area from the North pole to 60°S in latitude during the period of 1955-60 
was not so mnie different from those for 1950-55 and 1940-45. Generally speaking, 
the characteristics of the increase in 7 and decrease in X in the southern hemisphere 
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and the increase in X and rather decrease in Z in the northern hemisphere have been 
observed during these periods. These characteristics indicate that the decrease in 
moment of the earth’s magnetic dipole, which has been noticed for the past hundred 
years, was still continuing during the 1950-60 period. 

As the present isoporic charts for 1955-60 seem to be more reliable than those for 
1950-55, on account of the increase in the number of magnetic observatories, the 
present charts may be able to be compared with those for 1940-45 in certain details. 
By comparing these charts for the different periods, the westward drift of the foci of 
major isoporic vortices in each component can be noticed, the amount of the drift 
being ranged between 0° and 10° westwards over a period of 15 years. This problem 
will be discussed in a more analytical way in the following section. 

In the isoporic charts around the South pole, illustrated in Figs. 5, 7 and 9, re- 
markably large values of |Z! over and near the Antarctic may be noticed, amounting 
to about 2x10*y/year at their maximum. Generally speaking, in comparison with the 
North polar region, the geomagnetic secular variation is strikingly large and complicated 
in the South polar region not only for Z but also for X and Y, the annual rate of 
the variation in the latter region amounting to several times more than that in the 
former region. It has long been noticed that geomagnetic secular variations have 
been much more active in the Southern hemisphere than in the Northern hemisphere. 
Even so, it seems that the remarkable secular variation newly found in the Antarctic 
region is unexpectedly large, and it may need some particular interpretation to explain 
the existence of such large variations. 
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Fig. 1. Isoporic chart for X for 1955~60 in unit of 7/year. 
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Fig. 2. Isoporic chart for Y for 1955~60 in unit ANS 
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Fig. 3. Isoporic chart for Z for 1955~60 in unit oe r. Full circles show the positio 
of observing points and arrows represent Hack, Y). 


Geomagnetic Secular Variation during the Period from 1955 to 1960 91 


Fig. 4. Isoporic chart for X for 1955~60 in unit 
of 7/year in the Northern hemisphere. 
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Fig. 6. Isoporic chart for Y for 1955~60 in unit 
of ;/year in the Northern hemisphere. 
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Fig. 8. Isoporic chart for Z for 1955~60 in unit 

of y/year in the Northern hemisphere. Full cir- 

cles show the positions of observing points and 
arrows represent H=(X, Y). 
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Fig. 5. Isoporic chart for X for 1955~60 in unit 
of 7/year in the Southern hemisphere. 
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Fig. 9. Isoporic chart for Z for 1955~60 in unit 

of 7/year in the Southern hemisphere. Full cir- 

cles show the positions of observing points and 
arrows represent H =(X, ven 
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4. Spherical Harmonic Analysis of xX, Y and Z 


The isoporic charts given by Figs. 1-9 were subjected to the ordinary spherical 
harmonic analyses. In the practical procedures, readings were taken at points every 
15° latitude and every 30° longitude in cases where harmonic coefficients up to m=3 
are evaluated. These analyses were applied to xy Y and Z for their distribution 
within the area between 75°N and 75°S, and also for the area between 60°N and 60°S, 
the values at latitudes higher than 60° being ignored in the latter case. On the other 
hand, a more detailed analysis of Z, up to m=6, was carried out, based upon the 
isoporic charts between 80°N and 80°S, in order to obtain an underground equivalent 
current system in some detail representing the secular variation. In this case, readings 
were taken at points every 10° latitude and every 15° longitude. 

In all cases, results of analyses of the horizontal components, X and Y, were 
jointed together for determining the coefficients of the spherical harmonic series, while 
those for Z were treated independently of the horizontal components. Results of the 
above-mentioned analyses are given in Table II, together with those for 1950-55, 
obtained by Nagata and Rikitake, and for 1910-15, 1920-25, 1930-35, and 1940-45, obtained 
by Vestine et al, and in Table III. In these Tables, g,™ and h,™ are defined as 
usual in the following expression of potential variation ; i.e. 


: 


V=RYIDVe"cosmA+hrsinmd)P2(cosd) , (10% 


where FR denotes the earth’s radius. 

As will be seen in Table II, it seems likely that the general tendency of geomagne- 
tic secular variations represented by low degree and order harmonic coefficients has been 
presented nearly identical during the past half century. The main characteristics con- 
cerned are (a) decrease in the moment of the earth’s magnetic dipole, and (b) westward 
and northward shift of the position of the eccentric dipole. For the purpose of repre- 
senting these two characteristics more clearly, moment of the earth’s magnetic dipole 
(M), direction of the centred dipole (co-latitude 6) and longitude Ao of the North pole 
of the dipole field), position of the eccentric dipole (latitude ¢y and longitude 4, and 
distance r, from the centre of the earth in unit of the earth’s radius R) at the epoch 
of 1960 were computed, based on (A) 1955 values of the harmonics given by Finch and 
Leaton® and on (B) 1945 values given by Vestine et al. In the case of extrapolation 
from the 1955 values, the present results of secular variation for 1955-60 can be directly 
applied, while in the case of extrapolation from the 1945 value, the present secular 
variation data and those for 1950-55 are used. In Table IV, thus estimated values 
of M, 9, 40, mu, Ay and ry/R for cases (A), and (B) are given together with the cor- 
responding values at previous epochs. Now, it may be clear in the above results that 
M is still decreasing, and the eccentric dipole is drifting northwards, westwards and 
ae 5 a @o has been kept almost constant and the rate of decrease in 4) has been 


The speeds of the westward, northward and outerward drift (-A,, ur ru) of the 
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Table II. Coefficients of Spherical Harmonic Expansion of Geomagnetic Secular Variations 
for Various Epochs. (Unit=7/year) 
Coefficients | Author Vestine ef al Nagata-Rikitake Nagata-Syono 
1950-55 1955-60 

Bete: (Pe ee | ei ae 75° N-75°S 60°N-60°S 

X&Y Z OEE, omer 
g,° 25 28 23 9 9 6 15.0 16.0 14.6 1257 
git 1 4 1 2 0 2 12.0 6.4 11.4 6.6 
lore —7|/-7;,-—5 1 4 —1 | = 2.0 Les —4e 1 Meat 
go? — 7| —10} —14} -—18| —25 —=18. 1 22.3 | 25.3 | 90:7, || 220 
go! —1 1 1 0 0) SS | Cp) Sl Oak 7 le 187, 
h,! — 9} —14| —18| —20| —17 165) —6.0 217.5) | sore 720) 
go" 24 17 10 2 4 — 3 O200) = 0535) — 0529) — 024 
h,? —17 | —17 | —14] -—14] —14 => |) 16.8 || =17 295) =16.3' | =17.8 
g,° 6 4 4 2 6 —2 1,0} = 0.9 ih.) ORM 
2," —6)/—7|-—7 Fi 0 Ss OG Me ee a lee) 
h,! —6,/—57f—5)] = } 0 2) 7.4 8.4 5.9 7.8 
£," —1/—2 1 3 1 2 eS 0.7 ial 0.8 
bh, a ZS 4 4 7 —1 1-5 Aung) Te 4.1 
23° 13 7 3 0; —1 1 Aa ae 20 Sui 220 
bh? —11/ —11 | —13 | -—10} — 8 Sih 05S 58 


Table III. Coefficients of Spherical Harmonic Expansion of Z during 
1955—60. (Unit=7/year) 


[En™] 
PA: 0 1 2 3 4 5 6 
n 
| 13.6 5.9 
2 24.5 ie a).3 
3 1.9 i 1.1 —0.5 
4 — 4.4 rae | —2.0 =10 th) 
5 2.1 <0 0.2 1.0 23 =02 
6 eT, a 0.1 0.3 1.5 0.5 —0.6 
[h,™] 
ey 1 2 3 4 5 6 
n 
1 2.3 
2 —16.0 | gxdiz.i 
3 8.2 3.5 ah 7 
4 Oe. 120 3.5 eee 
5 1.8 7 a 0.7 0.2 
6 eye) S02 1.5 = vil 0.7 1.0 
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Table IV. 
Author of R 
Epoch analyzed data M 90 2o om Au ry/ 
emu 

1885 Schmidt‘? 8.346 x 1075 11°28" 29130: 5°36’N | 168°12’ 0.0461 
1922 Dyson-Furner“ 62169 11 34 290 54 9 13 161 16 0.0570 
1945 Vestine et al 8.082 ,, 11 24 289 54 13 56 154 05 0.0630 
1955 Finch-Leaton 8.068 ,, 11 42 291 03 15 41 150 49 0.0685 

Extrapolated from 
1960( A) 1. 1955 valties 8.046 ,, 11 30 290 34 16 40 149 20 0.0702 

Extrapolated from : 
1960( B ) Vootinels A948 values S020 4, 11 40 289 16 17 21 149 08 0.0679 


Table V. 
Average for 1955-60 Average for 1945-60 Average for 1922-55 
las +18'/yr +20'/yr +19'/yr 
Ou +12’/yr +14’/yr +12'/yr 
ry/R +3.4x10-4/yr +3.2x10-*/yr +3.5x10-4/yr 
M —4.4x10”"emu/yr —3.8x10"*emu/yr —3.1x10"emu/yr 


eccentric dipole and the annual rate of change in the dipole moment (M) are estimated as 
given in Table V. So far as the rate of secular change in moment M and that of the drift 
of the eccentric dipole during recent years are concerned, it seems likely that those values 
estimated directly from the secular variation data obtained at a number of geomagnetic 
observatories may be more reliable than those estimated from the difference between 
magnetic charts at different epochs, which are compiled by different authors with diffe- 
rent methods. 

As for the drift of the geomagnetic non-dipole field along latitude circles, each 


—0.2 0 0.2 0.4 0.6 
_ or degree|yr. 
ot 


Fig. 10 Histogram of drift velocity of spheri- 
cal harmonics (254) of the geomagnetic 
secular variation along latitude circles. 
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spherical harmonic term (g”, h”) of secular variation distribution can be examined 
directly. If the non-dipole field potential represented by 


W=R3} S\(gucosma+hrtsinma)Px(cosd) (2%) 


n=2 7 


has been drifting along latitude circles without changing its distribution form, its change 
is given by 
_ OV AA 


~ Oh “7 RE Lmlhrcosmi— gusinma]e 7 Pn(cosd) , (3) 


where 0A/0¢ gives the velocity of the westward drift. Then, from (1) and (3), we get 


Westward-drift velocity=“*=-&" —— hin : (Cb) 
ot mh™ mg, 
From the values of (g”, h™ given in Table III, and the values of (g”, h™) for 1950 
obtained by Finch and Leaton, numerical values of 04/0t were computed for 18 couples 
of coefficients for 4=n=>2. The results are shown in Fig. 10, as a histogram of 024/0t 
values. It may be seen in this diagram that each of the harmonics is drifting mostly 
westwards with a velocity of several decimals of a degree per year. 


5. Equivalent Current System of Geomagnetic Secular Variation over the 
Surface of the Earth’s Core. 


From the harmonic coefficients, gn and h”, for n<6, given in Table III, the equiva- 
lent current system on the surface of the earth’s core, namely on the spherical surface 
of 2900 km in depth from the earth’s surface, can be obtained. Figs. 11-13 illustrate 
the distribution of the equivalent current system, where the harmonic coefficients “=7 
are ignored. 

It seems that the distribution of secular variatent can well be represented by the 
equivalent current system of <6 over almost all regions over the earth except the 
Antarctic region. As already mentioned, regional anomalous variation around the south 
pole is extremely large and local, so that a difference (observed value minus computed 
value by this equivalent system) still amounts to from —40y/yr to 407/yr for Z in this 
Antarctic region. Attention may be especially paid to a pair of Z foci in the vicinity 
of East Antarctica, a positive Z focus on inland in East Antarctica and a negative Z 
focus near Heard Island, the distance between them along the great circle and the 
difference between their Z foci is the result of a pair of radial magnetic dipoles on the 
earth core’s surface, the position and intensity of these dipoles are estimated as illustrated 
in Fig. 14, which shows a cross section of the earth along the great circle passing 


through these foci. 
As shown in the figure, the upward magnetic dipole has about 3.5 10”emu/yr in 
its annual rate of change in moment, while the downward one about 15° apart is about 


2/3 in its change of moment. The two curves of the upper half of the figure illustrate 


the distribution of Z and H derived from the pair of dipoles and the circles represent 
The above-mentioned result 


the corresponding observed values at the four stations. 
f a pair of vortices, one 


may indicate that a change in electric currents of the form o 
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Fig. 11. Equivalent current system of the geomagnetic ait variation at depth of "2900 km 
in unit of 10 Ampere/year. 


Fig. 12. Equivalent current system of the geo- Fig. 13. Equivalent current system of the geo- 
magnetic secular variation in the Northern magnetic secular variation in the Southern hemi- 
hemisphere at depth of 2900 km in unit sphere at depth of 2900km in unit 
of 10 Ampere/year. of 10* Ampere/year. 


clockwise and another counter-clockwise, is taking place near the core’s surface beneath 
this region. This means that equivalent currents between the two foci of vortices on 
the core’s surface are changing at the rate of about 0.1 Amp/cm/yr. (In Fig. 13, con- 
cerned only with lower spherical harmonics of n<6, the rate of change in electric current 
density amounts to about 0.05 Amp/cm/yr.) It may be said that this rate of change 
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Fig. 14. Distribution of Z and H along the 
great circle plane passing through Syowa and 
Heard Island in the Antarctica. | 
Above: observed values of Z and H and 
those produced by a pair of magnetic 
dipoles illustrated below. 

in electric currents is significant, compared with the equivalent current of the earth’s 
main dipole field. 

Intensity of electric currents over the core’s surface, equivalent to the earth’s dipole 
field, is approximately given by J,=1.4cosy,, Ampere/cm, where ¢,, denotes geomagnetic 
latitude, and therefore their intensity in this region-where ¢,~60~70° is about 0.5 
Ampere/cm or less. Then, the total amount of current change of the vortex-form in 
the Antarctic area over several years ought to become as extensive as the intensity of 
the dipole-field current. On the other hand, observed evidence has shown that the large 
rate of geomagnetic secular variation in this area has continued throughout the past 
30 years, at least®. 

It seems unlikely therefore that the intense regional secular variation of a pair of 
vortex represents in form the simple fluctuations in the electric current system within 
the core responsible for the geomagnetic dipole field owing to its instability. Alternati- 
vely, it might be rather reasonable to presume that the change is caused by the break- 
ing out of a part of the inner toroidal magnetic field, such as is assumed by Bullard, 
or by the flowing up and breaking out of a ring of the toroidal magnetic field near the 
core’s surface, such as is assumed by Alfvén™ in explaining a pair magnetic field of 


the sunspots. 
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Abstracts of the Papers Presented at the 28 th Annual Meeting, 
Kyoto, Oct. 30-Nov. 1, 1960. 


OGAWA T. and SaGa §S._ (Geophysical 
Institute, Kyoto University) Charge on Rain 
Drop (Continued report )——The problem on 
the electricity of precipitation is to search the 
electric charge on rain drop for the history of 
charging. This paper is our primary attempt 
to know the magnitude of the original charge 
on rain drop by considering the rain current 
at the ground on the standpoint of Wilson’s 
ion-capture theory. In general, the rain cur- 
rent is represented by the equation 
i= | NQVax, where N is the size distribution 
of rain drops, Q@ the charge and V the falling 
velocity of a rain drop. All of these variables 
are the functions of drop size x. If the charge 
on rain drop can be negligible small compared 
with 3Xa? (X: the electric field in esu, a: the 
radius of rain drop in cm) at any instant dur- 
ing the time interval in which the rain drop 
falls from a cloud to the ground, the charging 
rate of rain drop is represented by dQ/ct= 
—3z1Xa* according to the calculation of Whipple 
and Chalmers, where Z is the polar conductivity 
in esu. Applying the result of integration of 
this equation to the rain current, we tried to 
dissect Simpson’s empirical formula. As the 
result of this dissection, a relation of n~695/ 
is obtained, assuming the original charge on 
rain drop Q,=3nXa*. The relation is given in 
the following table for some values of 2%. 


Z (esu) ) n | 
; 
1.5 10-5 (Average value at Kew) 0.0104 | 


110-4 (Average value at land | 
stations) 0.0695 | 


2x 10-4 | 0.139 


Thus rain drop starts a cloud with a small 
charge compared with 3Xa’ in the same sign 
as the electric field, and reverses its sign during 
the fall to the ground. As m must be a certain 
definite value, the process in which the charge 
on rain drop changes its sign in the region 
below cloud proceeds depending on the value 
of conductivity; this may be the reason why 
the mirror image effect has been found clearly 
at some stations and not so at others. 


SAGA S. and OGAWA T. (Geophysical 
Institute, Kyoto University) Electric Field 
below Cloud On the disussion of the charge 
on rain drop, there remains a serious problem 
on the electric field below the cloud, that is the 
discrepancy between the Altielectrograph results 
by Simpson and Scrase and Simpson and 
Robinson, and theoretical calculation by 
Whipple and Scrase on the basis of Wilson’s 
Space Charge Theory. Taking the life time 
of small ions into the consideration, the 
authors found the possibility of rapid in- 
crease of the field below the level of 100m; 
the magnitude at the level of 100 m is about ten 
times as large as that at the ground. This 
rapid increase may have caused the erroneous 
estimation of the coefficient in the relation be- 
tween the current and the field strength used by 
Simpson and Scrase, and so the above discre- 
pancy may have arisen. Standing on this sup- 
position, the small ions produced by point dis- 
charge do not ascent above the height of a few 
hundred meter, and the electric field remains 
nearly constant above that level. Consequently 
it may be thought on Wilson’s process below 
the cloud; that the rain drop captures hardly 
ions produced by point discharge, but mainly 
small ions which contribute to the normal con- 
duction current under the high electric field. 


MISsSAKI M. (Meteorological Research In- 


stitute) Studies on the Atmospheric Ion 
Spectrum II In the previous investigation 
on the mobility spectrum of atmospheric ions 
with mobilities greater than 0.4 cm?/volt sec, it 
was deduced that small ion concentration was 
too small to account for the electrical conduc- 
tivity, so that we were forced to conclude that 
there might be appreciable contribution of large 
ions to conductivity. This is likely to occure 
when the air is polluted as in the case of the 
large city. To leave little room for doubt about 
this, however, we must further proceed with 
the studies on the spectrum in the region smaller 
than 0.4 cm*/volt sec. A thick probe, 12cm in 
diameter, was produced for the measurement of 
spectrum in the mobility region down to 0.002 
cm?/volt sec. This probe has the same length 


(9 ) 
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with the slender one for small ions, and can 
be replaced with it into the aspiration cylinder. 
In the present works, spectra were obtained 
between 0.4 and 0.02 cm?2/volt sec. Analysis of 
several examples showed that the contribution 
of large and intermediate ions to the conduc- 
tivity attains to a few tens percent of total con- 
ductivity. This is in agreement with the pos- 
tulate of previous works. 


UCHIKAWA K. (Japan Meteorological 
Agency) On the Atmospheric Electric 
Phenomena in the Upper Air (Variations 
in the Vicinity of Exchange Layer) 
Measurements of electrical conductivity and 
potential gradient in the upper air using special 
designed radiosondes were carried out at four 
stations in Japan (Sapporo, Tateno, Hachi- 
jojima and Kagoshima) during IGY and IGC. 
An investigation of these results shows that 
height of the exchange layer has its mean value 
of 2200 m, and varies slightly from season to 
season, showing high in summer and low in 
winter. Considerable discontinuities of columner 
resistance and space charge are observed at the 
boundary between the exchange layer and the 
upper level. This discontinuity, however, dis- 
appears when meteorological disturbances (cy- 
clone, front etc) pass through the observation 
point, and then continuous distributions of these 
elements last for several days. Calculating from 
observed values of temperature and wind at the 
level of the boundary, descending current (or 
ascending current) is presumed when the ex- 
change layer exist (or not). 


TAKAGI M., (The Research Institute of At- 
mospherics, Nagoya University) On the Dis- 
charges in the Thundercloud. The dis- 
charge in the thundercloud is composed of a 
main slow discharge process and of many rapid 
local breakdowns occurring in a slow discharge 
process. The investigation of the polarities of 
electro-static field changes and of the amplitude 
relations between static field changes and light 
emissions, both produced by such local break- 
downs, has led to the following results. The 
cloud discharge is generally produced between 
the main upper positive and lower negative 
charged regions in the thundercloud, and the 
successive positions of local breakdowns are 
descended with the downward development of 
a slow main discharge process started from the 
positive charged region, but a local breakdown 
is accomplished in itself by the negative short 
streamer advancing rapidly upward. While the 


discharge inside the cloud, appearing between 
and after the multiple ground strokes, is similar- 
ly constructed by a slow positive discharge pro- 
cess and rapid negative local streamers, except- 
ing that the directions of advances are respec- 
tively contrary with those of the cloud discharge. 


TAKEUTI T. and IsIKAWA H. (The 
Research Institute of Atmospherics, Nagoya Uni- 
versity) Ground Discharge Mechanism 
Deduced from Electric Field Records Ob- 
tained Simultaneously at Two Places. 
From electric field records obtained at Maebashi 
and Shibukawa, which is north apart 13 km 
from the former, two results about ground 
discharge are obtained. (1) We obtained a few 
couples of electric field change records by cloud 
discharges preceding the first ground strokes, 
which is named pre-preliminary discharge in 
our paper.* As the number of records is too 
poor to discuss about it, we only show those 
records. (2) It has reported that ground strokes 
neutralize successively upward a _ negative 
charged column in a thunder cloud. We have 
deduced that the negative charged column in- 
clines, so that the stroke pass in the cloud has 
horizontal direction component, and the J and 
F processes named by Malan also have horizontal 
direction component. *Takeuti T. etc. Proc. 
Res. Inst. Atmo. vol. 7 p. 1 1960. 


ISIKAWA H. (The Research Institute of 
Atmospherics, Nagoya University) The Inter- 
pretation of the Pulsive Structure of the 
Electromagnetic Waveform of a Prelimi- 
nary Discharge——The time intervals between 
radiation pulses appearing on the waveform of 
a preliminary discharge have been measured on 
electromagnetic waveforms. The variation of 
time intervals between appreciable pulses with 
the progress of a preliminary discharge mea- 
sured on the waveforms fairly well coincides 
with the step interval variation with the pro- 
gress of a stepped leader measured on the photo- 
graphs of stepped leaders. Numerous small 
pulsations superposed on a series of appreciable 
amplitude pulses, however, usually have the 
frequency of repetition more than 50 kc and 
can not be considered to correspond each ‘to a 
step streamer. The variations of pulse intervals 
with the discharge progression have been com- 
pared between a type and f type stepped leaders, 
however, we can not find any essential differ- 
ence existing between the two. The difficulty 
in photographing a clear f type stepped leader 
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therefore should be attributed to the decrease 
in the streamer luminosities. 


a URITZT OO: (Electrical Communication 
Laboratory, Nippon Telegraph and Telephone 
Public Corporation). Unified Theory for 
Tropospheric Propagation of Radio Waves 
——A theory of radio propagation in a homo- 
geneous atmosphere valid for all regions within, 
near and beyond the radio horizon is presented. 
Use of the saddlepoint method which proves 
useful for analysis of the normal propagation 
within the horizon reduces the problem to find- 
ing an appropriate root of the saddle-point equa- 
tion. It is found that the magnitude of this root 
is closely related to the field strength received 
at a distance from the transmitter and that its 
imagi nary part, in particular, determines the 
nature of the received field strength which varies 
markedly from a region to another. Depen- 
dence of the received field strength on the 
gradient of atmospheric refractive index near 
the earth surface is discussed. 


Miya K. 
Co., Tokyo) Bearings of Atmospherics at 
HF Bearings of atmospherics at 5, 10 and 
15 Mc/s were measured at Ohira near Tokyo 
in December, 1955, by using a direct vision 
type direction finder. Ranges of bearings were 
obtained from a photograph taken with 1 minute 
exposure. Atmospherics generally arrive from a 
wide range of southerly bearings except a range 
between 320° and 50°. Atmospherics at 5 and 
10 Mc/s change their bearings clockwise, i.e. 
from south to northwest in the morning and 
from northeast to south in the afternoon. On 
the other hand, atmospherics at 15 Mc/s almost 
arrive from southerly bearings throughout a 
day. The above characteristics of bearings are 
explained by differences of strengths of atmos- 
pherics induced in beam antennas pointing to 
different directions of 53°, 244° and 337°. It 
is shown that these results can be well ex- 
plained by using the worldwide mop of atmos- 
pherics, which is rearranged from the maps of 
CCIR report, in the expression of U.T. 


(Kokusai Denshin Denwa 


TsSUBOKAWA I., FuyirA N. and SETO 
T. (Geographical Survey Institute) On Proton 
Precession Magnetometer (Comparison 
with G.S.I. Magnetometer)—— The com- 


ponent proton precession magnetometers report- 
ed by Tsubokawa I. in 1957, are now being 
operated in the Kandzan Geodetic Observatory, 
recording H and F continueously (recording 
speed: 50 mm/h, sensitivity: about 1 mm/y, time 
interval: 30 sec.). The head of the magneto- 
meter for H is placed in the field of Fanselau 
coil whose current is regulated in the accuracy 
of 0.5 % automatically. The coil axis error and 
the level error are being checked and corrected 
once a week. The dip computed from the 
records of F and H of proton precession mag- 
netometer are compared with that obtained 
directly by G.S.I. magnetometers No. 27 and 
No. 41 which are electromagnetic magneto- 
meters with rotating coil detector. These in- 
struments are respectively placed in different 
stations. The difference of the geomagnetic 
field among the stations where the proton pre- 
cession magnetometers and G.S.J. magneto- 
meters are placed, are determined by above 
mentioned G.S.J. magnetometers, exchanging 
them alternately. The result is as follows: 
Ip—Ig=0.'00-+40.'03 

Where Ip is the dip computed by the proton 
precession magnetometers and Ig is the dip 
obtained by G.S.I. magnetometers. This shows 
Ip and Ig agree with each other within the 
accuracy of observation. 


Kato Y,. and Other Members of Research 
Group of Air borne Magnetometer (Geophysical 
Institute, Tohoku University) On the Air- 
Borne Magnetometer (For the Measure- 
ment of the Vertical Component We 
constructed the new type air-borne magneto- 
meter for the direct measurement of the 
vertical component of the earth magnetic 
field. The detector core of this magneto- 
meter is fixed on a_ vertical air-gyro- 
scope which holds the vertical direction by 
means of high speed rotation. The method of 
measurement of the vertical component of the 
earth’s magnetic field is an application of the 
Tohoku University Method for the measure- 
ment of the week magnetic field and one kind 
of so-called saturable core magnetometer. 
The magnetic core of easily saturable ferro- 
magnetic material of high permeability is set 
on the vertical axis of revolving body of the 
gyroscope, and the primary and secondary coils 
are set on the bracket of the gimbal (on which 
the gyroscope is mounted). The rotating body 
of gyroscope having a large angular momentum 
does not make any change evenwhen the gimbal 
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receive some acceleration in a short time. Thus 
the axis of the gyroscope or the magnetic core 
of the magnetometer hold the vertical direc- 
tion, and we can detect the vertical component 
of the earth magnetic field in the airplane. 
We could observed the magnetic anomaly 
above 3000 m of Mt. Mihara, Ooshima and we 
could separate the magnetic anomaly due to 
the subterranean structure of the Mt. Mihara 
and the basement structure of Ooshima. 


HARADA Y. (Geographical Survey In- 


stitute) Geomagnetic Anomaly Distribu- 
tion in Japan From the results of the 
second order magnetic survey at about 800 
stations established during the periods from 
1952 to 1959 in Japan, magnetic charts of 
three components for 1955.0 are obtained. 
Method of the least square with a _ condi- 
tion equation, rot H=O, is used for the seven 
sectorial zones independently, and then magnetic 
chart is combined from the seven partial charts. 
A) From the main distribution represented 
by second order polynomial equations, two 
anomalous zones for Y-component are found 
along Fuji volcanic range in Honsyu and Hidaka 
Mt. range in Hokkaido. The magnitudes of 
the anomalies are estimated to about 20’ in 
declination. B) From the residual charts, (a) posi- 
tive anomalies of 4z are found in alluvium and 
pleistocene districts; (b) negative anomalies of 
4z are found in the mountaineous districts, 
Chubu, Shikoku, Chugoku and Kyushu; (c) 
large positive anomalies in 4z over than 2007 
along the volcanic ranges are remarkably found; 
(d) generally speaking, Western Japan is not so 
much disturbed, while Eastern Japan is much 
disturbed. 


YUKUTAKE T. (Earthquake Research 


Institute, Tokyo University) Preliminary 


Study on the Secular Variation of the 
Geomagnetic Field Time variation of the 
magnetic field at the earth’s surface consists of 
two different kinds of variation ; one originated 
in the steady movement of the field and other 
caused by the growth or decay of the drifting 
field. Some trials to evaluate their contribu- 
tions to the field variation have been com- 


menced. Spectrum analysis of the varying magne- 


tic field based on palaeomagnetic investigations 
showed that phenomena having such periods as 
7000, 3000, 1800, 1200, 800 and 400 years were 


conspicuous. The most of these rather seem 
to be associated with the drift of the main 
features of the field at the earth’s surface than 
the variation related to the vissicitude of the 
fluidal motion in the core. Analysis for short 
intervals of recent observations revealed the 50- 
and 25-year period, though whether these 
originate within the earth or not is unknown 
yet. Assuming that the centered dipole has 
shown no appreciable movement over thousand 
years, the variations in declination and inclina- 
tion caused by the steady westward drift of the 
non-dipole field were constructed for last one 
thousand years in Japan and compared with 
archaeomagnetic results. Good correspondence 
can be seen regarding the main features of 
variation when drift velocity of about 0.35°/ 
year is assumed. The same kind of comparisons 
suggest that the drift velocity in the past might 
be faster than now. 


NAGATA T. and SYONO Y. (Geophysical 
Institute, University of Tokyo) Geomagnetic 


Secular Variation during 1955-60 Geo- 
magnetic isoporic charts of X, Y and Z for the 
period of 1955-60 were constructed, based upon 
observed data at 99 stations including 12 An- 
tarctic ones. The results were subjected to 
spherical harmonic analysis. Main points of 
the results are as follows. 
(a) Moment of the earth’s magnetic dipole 
has continued to decrease. The values of the 
moment () and position (@, 4) of north pole 
of the dipole at 1960.0 are :— 
M (1960)=8.046 x 10> emu, 
%o(1960)=78°30'N, 2)(1960)=290°34’. 
(b) The geomagnetic eccentric dipole has 
continued to drift westwards, northwards and 
outwards. The position (Lat.=yy, Long.=2Zy, 
distance 7, from the earth’s centre in unit of 
the earth’s radius R) of the eccentric dipole at 
1960.0 are :— 
gu=16°40’N, 2% =149°20’, r/R=0.0702, 
while its velocity of drift is given by 
ou=+12"/yr, 2y=—18'/yr, 
ry=+3.0x10-*R/yr. 
(c) Velocity (A) of the westward drift was 
estimated by 
A=—g/mh and A=+h/mg for harmonics 
of- mSnS3. Average of A values is 
—15’/yr. 
(d) There is a twin of extremely intense 
positive and negative foci of Z near the coast 
of East Antarctic. There intensity amounts to 
about 200 7/yr. 
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KONDO I., MuRAYAMA T., Mori S., 
OKADA H., MAKINO T., SAKAKIBARA S., 
KATO S and Honzawa T. (Physical Institute, 
Nagoya University) On the Observation of 
High Energy Cosmic Ray Particles by 
Cosmic-Ray Telescope No. 3 (I) Detec- 
tion Efficiency Calculations were made 
on the detection efficiency of cosmic ray 
particles for newly built Cosmic-Ray Telescope 
No. 3 (J. G.G. 12, 46, 1960). (1) Maximum 
number of photo-electrons available at the 
photo-cathode of multiplier photo-tubes was 
determined as 34 for each passage of a 
particle, including the effect of the wave 
length shifter (7;=34). (2) To suppress 
accidental coincidences to several counts per 
hour, photo-tube pulses are discriminated at 
such a pulse height that corresponds to 1.5 
electrons at the photo-cathode (m,=1.5). 
(3) With above mentioned values of 2, and 
m,, the solid angle for each four-fold coin- 
cidence can be estimated to be 4°x6°, and 
about 7° in diameter for a ‘‘ Y-element’’ (mix- 
ture of three 4~f coincidences). The total field 
of view is about 17° in diameter which consists 
of six ‘‘ Y-elements’’ and six ‘‘ V-elements’’ 
(mixture of two 4-/ coincidences). As these 
elements overlap one another, incoming particles 
are detected with high efficiency (>95%) in 
this field view. (4) The effective area of 
the telescope is about 10m?. (5) Expected 
intensity of vertical cosmic rays to be 
observed by this telescope is estimated 
with the knowledge on the energy spec- 
trum of vertical »-mesons and electrons, amount 
of Cerenkov light as a function of a particle 
velocity and the geometry of the telescope 
(shown in (3) and (4)). The deduced value 
(4,000 per hour for each ‘ Y-element”) is in 
good agreement with the observed one (3,600 
per hour). 


Konpbo I., MURAYAMA T., Mori S., 
OKUDA H., SAKAKIBARAS., MAKINO T. and 
HoNnNzZAWA T._ (Physical Institute, Nagoya 
University) On the Observation of High 
Energy Cosmic Ray Particles by Cosmic- 
Ray Telescope No. 3. (II) Zenith Angle 
Dependence Using the newly built cosmic- 
ray telescope, directional intensity of high 
energy cosmic-ray particles (mu-mesons Se lO 
GeV, electrons > 150 MeV) was measured. 


In the course of measurement of zenith 
angle dependence, two sources of background 
pulses were noticed. J.e., (1) pulses due to 
Cerenkov radiation in photomultiplier windows 
emitted by local shower particles. (2) Accidental 
coincidence pulses by noise of photomultipliers. 
(1) is much more important (~1/min, for a Y- 
element) and will be reduced by anticoincidence 
method while (2) has smaller effect (~0.1/min., 
N=2.10'/sec, c=1.5. 10-7 sec). The counting 
rate of the telescope was corrected for 
these background effects and the zenith angle 
dependence of the intensity was obtained. 
This was in good agreement with theoretically 
estimated zenith angle dependence of high 
energy cosmic-ray intensity. From this zenith 
angle dependence, it was found that 50 
percent of primary cosmic-ray particles which 
produce these high energy cosmic-ray particles 
have energies between 30 and 80 GeV for z=0° 
and 170 and 330 GeV for z=80°. The east 
west asymmetry of cosmic-ray intensity was 
found to be less than 3 percent for the energy 
range observed by this telescope, even at 
large zenith angle. Some preliminary results 
of continuous observation of cosmic-ray in- 
tensity by this telescope was also presented. 


KAMIYA Y., SAGISAKA S., UENO H. and 
KATO S. (Nagoya University) On the Con- 
struction of Cosmic-Ray Telescope with a 
Solid Iron Magnet In order to determine 
the sign of charge and the momentum in high 
energy cosmic-ray particle, a mew cosmic-ray 
telescope with a solid iron magnet (Cosmic-Ray 
Telescope No. 4) is now under construction at 
Nagoya University. The design of this appara- 
tus were reported in the preceeding paper pre- 
sented the annual meeting, May 16-18, 1960 
(published in J.G.G. Vol. XII, No. 1). This is 
preliminary report. The total magnetic flux in 
this magnet was measured by means of search 
coil and fluxmeter. The mean induction is 20 
kilo gauss with excitation current 30 ampere 
and the variation of the excitation current is 
about +1%. The efficiency of the neon hodo- 
scope is about 80%. 


KITAMURA 7T. (Geophysical Institute, 
Kyoto University) and KODAMA M. 
(Institutd of Physical and Chemical Research) 
Time Variations of Cosmic Ray Intensity 
in the Antarctic Region (I) Atmospheric 
Temperature Effect and Diurnal Variation 
——Cosmic-ray observations at Syowa Base, 
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Antarctica, have been carried out by using a 
Neher type ionization chamber during 1957-58 
and a cubical telescope having two plastic scin- 
tillators since 1959, respectively. The amplitudes 
of diurnal variation of cosmic-ray meson com- 
ponent obtained from these observations are 
calculated in every month together with those 
at Mawson and Tokyo. It is seen from the 
figure (not given here) that the month-to-month 
variations of the amplitudes at the two stations 
in the Antarctic region are obviously different 
from that in the middle latitude station. This 
discrepancy seems to be due to the difference 
of the atmospheric temperature variations bet- 
ween the higher latitude and the middle. 
In order to extract the parts due to the 
atmospheric temperature effect only from the 
curves in Fig. 1, it was firstly corrected for 
the primary cosmic-ray variation by using the 
data of neutron monitor at Mawson. Then, 
the temperature coefficient of meson component 
was deduced from the several correlations bet- 
ween the different heights of isobar level and 
cosmic-ray intensities. In the results, the best 
correlation coefficient was obtained not in 300 
mb, above which the best correlation in the 
middle latitude is usually found, but in 50 mb. 
Using the temperature coefficient thus ob- 
tained, 5%/km, it is possible to presume the 
amplitudes of diurnal variations of 50 mb isobar 
level, that is, they are about 12 meters in 
summer and about 7 meters in winter. 
To ascertain the above mentioned fact from 
the observational data of Radio Sonde, the 
diurnal variations of isobar levels at Tateno 
(Japan) and Kap Tobin (Greenland) are illust- 
rated. There are no apparent seasonal changes 
in both 100 mb and 300 mb at the former sta- 
tion, whereas such a seasonal change is obvious 
at the latter corresponding to the same latitude 
as Syowa Base. Therefore, it is reasonable to 
say that the apparent seasonal change of diur- 
nal variation on cosmic-rays at higher latitude 
are mainly due to the atmospheric temperature 
effect. 


KITAMURA M. and Tezuka T. (Me- 
teorological Research Institute) Qn the Daily 


Variation of the Cosmic-Ray Intensities 
(II) Data of cosmic-ray intensities observed 
by the neutron-monitors distributed all over 
the world during IGY are analysed to obtain 
the latitude dependence of diurnal and semi- 
diurnal variations of the cosmic-ray intensities 
for all days during IGY and storm times, 


respectively. The results are as _ follows: 
(1) The amplitude of the diurnal variation of 
cosmic-rays in the auroral zone for both all 
days and storm times are larger than those of 
other regions. Furthermore, increments of the 
amplitudes of the diurnal component are larger 
at high latitude than other latitudes. (2) The 
amplitude of the semi-diurnal variation shows 
the maximum value in the equatorial zone for 
both all days and storm times. And the energy 
dependence of increments of the amplitudes of 
the semi-diurnal variation during storm times 
are not so remarkable. (3) The maximum 
time of the diurnal variation for all days is about 
Ilh at the equatorial zone and later in the high 
latitude. They are shifted towards the early 
morning in the low and middle latitude during 
storm times and are delayed in the auroral 
zone. (4) The maximum time of the semi- 
diurnal variation is about Oh at the equatorial 
zone and later in the high latitude. They are 
delayed in the high latitude during storm times. 


KANNO T. (Department of Physics, Fuku- 
shima University) On the World-Wide Dis- 
tribution of the Daily Variation of the 


Cosmic Ray Neutron Intensity and its 
Variation (II). Tm, the local time of the 
maximum intensity in the diurnal variation of 
the cosmic-ray neutron intensity, was derived 
from the data of 39-43 stations during the 
period of I.G.Y. It seems that there are the 
latitude effect and the longitude effect in ary 
Tm in the vicinity of the geomagnetic equator 
advances and the magnitude of the latitude 
effect becomes large as the solar activity be- 
comes intense. The longitude effect varies 
every month and the magnitude of its variation 
distributes with a certain tendency in the world. 
It may be suggested that the former is due to 
the modulation with magnetic clouds in the 
interplanetary space, and the latter is due to 
the deformation of the lines of force of the 
terrestrial magnetic field by the solar magnetic 
clouds. 


OsHIO M. (Faculty of Science, Tokyo 
University) KODAMA M. (The Institute of 
Physical and Chemical Research) World-Wide 
Distribution of the Cosmic-Ray Neutron 
Intensity and its Time Variation during 
Magnetic Storms In order to study physi- 
cal state of magnetic cloud during magnetic 
storm, using the bi-hourly values of cosmic ray 
neutron intensities obtained at 35 stations during 
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IGY, intensity variations in the directions fixed 
to 12 local times are analyzed, where each 
station is divided into four groups according to 
the geomagnetic cut-off rigidity. At first, March 
25, 1958 cosmic ray storm is adopted as an 
example not to be accompanied with large 
decrease in horizontal intensity of magnetic 
storm. Analysis seems to show that the ampli- 
tude of anisotropy of primary cosmic-ray 
particles reached the earth varies damping 
with 3~4 day period, and its direction becomes 
opposite in the latter half of the storm. 


TORIZUKA K. (Tokyo Gakugei University) 
and WapaA M. The Institute of Physical 
and Chemical _ Research, Cosmic-Ray 


Storms in July 1959——Cosmic-ray storms 
occurred on 15 and 17 July 1959 were studied. 
Those variations were recorded by_means of a 
high sensitive meson detector installed at the In- 
stitute of Nuclear Study, Tokyo. The observation 
error was 0.1% of the mean intensity for the total 
count of several minutes. The same kind of 
observation was also carried out at Massatusetts 
Institute of Technology, U.S.A... A comparison 
of the data of the two stations resulted the 
followings. (1) The variations of the 15 event 
at the two stations were in good parallelism, 
whereas those of the 17 event were in large 
difference with each other. (2) Two well sepa- 
rated increases were found during the decreas- 
ing period at 13-16 hr U.T., 15 July simulta- 
neously at the both stations. The variations 
were: The first minimum was 1333 U.T., the 
first maximum was 1357, the second minimum 
was 1445, and the second maximum was 1506 
in which the error of each time was 6 minutes. 
(3) The difference of the 17 event may be 
interpreted partly by the existence of many 
short lived increases during the period when 
the intensity was depressed. The increases 
were observed only at MIT which caused appre- 
ciable anisotropy of the intensity of coming 
particles. 


KODAMA M. (The Institute of Physi- 
cal and Chemical Research) A Large Fluc- 
tuation of Cosmic Ray Intensity on 
July 18-19, 1959——Cosmic ray intensity 
largely fluctuates during the recovering 
period 18th-19th of Forbush decrease on July 
17, 1959. Cosmic ray stations represent such 
a fluctuation are not all over the world but 
limited to the geomagnetic latitude range 40°- 
70° and longitude range 20°-200°. This phenom- 


enon occurs in not only the northern hemi- 
sphere but also the geomagnetic conjugate region 
in the southern hemisphere. It seems that Van 
Allen belt has a close relation with the present 
event. 


TAKAHASHI H. (Faculty of Liberal Arts, 
Iwate University) Latitude Effect of Cosmic- 
Ray Storms Associated with Magnetic 
Storms——tThe statistical result about depen- 
dency of the latitude effect of cosmic-ray storms 
associated with magnetic storms on the magnet- 
ic storms is stated. The world-wide cosmic- 
ray data (14 stations) on mountain altitude and 
geomagnetic data obtained at Kakioka Magnetic 
Observatory during IGY are used. The result 
obtained is that the latitude effect of cosmic- 
ray storms associated with magnetic storms 
depends on the magnetic storms and it is great- 
er on the large magnetic storm than on the 
small one. The further study will give the 
detailed discussions. 


Kupo S. and MuRAKAMI K. _ (Inst. 
Phys. Chem. Res.) On the Relation bet- 
ween the Cosmic-Ray Diurnal Variation 
and the Other Phenomena ——It is very 
important to find out the substantial phenom- 
enon related with the diurnal variation of the 
cosmic-ray intensity. In the present report, the 
frequency distributions of the amplitude of the 
diurnal variation, associated with the geo- 
magnetic or solar events, were compared each 
other on the shape of the distribution. After 
the comparison of the distribution during a 
certain period, it is concluded that Type IV 
outburst in solar flare is a very essential phenom- 
enon for the classification of the diurnal 
variation, rather than the cosmic ray storm or 
the other phenomena. And the change of the 
shape of the distribution was convinced with 
the lapse during one and half years. This 
change also is explainable, only if consider the 
change of the distribution of the diurnal varia- 
tion associated with Type IV outburst. 


TAKAHASHI H. (Faculty of Liberal Arts, 
Iwate University) Correlation between Time 
Intervals from Onset of Cosmic-ray Storm 
to the Least Intensity and Solar Activity 
_—There is rather good correlativity between 
the yearly mean (t;) of time intervals from 
onset of cosmic-ray storm to the least intensity 
and the yearly mean (R) of relative sunspot- 
numbers. f; increase with increasing R, and 
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they indicate a variation with solar cycle. This 
fact is interpreted by a simple solar cloud model. 


Kamtiya Y. (Nagoya University) Cosmic- 


Ray Storms and Type IV Solar Radio 
Outbursts In the preceeding paper (Rep. 
of Iono. and Space Res. in Japan, vol. XIII, No. 
2, 1959), relations between the type of solar 
radio outburst, cosmic-ray storm and magnetic 
storm were studied using the data obtained 
during IGY, and the following results were 
obtained. (1) Almost all cosmic-ray storms were 
associated with Type IV solar radio outbursts. 
(2) The size of a cosmic-ray storm is independent 
of the meridian distance of the radio outburst, 
but that of a magnetic storm depends on it. 
In this paper, these results were more con- 
firmed by using these data of the years 1959 
and 1960. 


SAKURAI K. (Geophysical Institute, Kyoto 


University) Solar Protons and Electrons 
—We examine the behavior in the solar outer 
atmosphere and the geomagnetic dipole field of 
high energy protons and electrons which are 
produced from the sun: (a) we calculate, in 
detail, the motions of the solar cosmic-ray pro- 
tons with a few Mev—several 100 Mev energies 
in the geomagnetic dipole field and then, seek 
for their precipitation lines on the earth’s sur- 
face. In consequence, we obtain the following 
results; (1) there is anisotropy on their in- 
cidence, (2) their precipitation lines show the 
circular ones around the geomagnetic north 
pole, (3) their incident region is restricted in 
the polar cap region only, and (4) statistically, 
their incidence tends to concentrate to the morn- 
ing hemisphere. Since the observed facts 
coincide well with the theoretical calculation, 
we may consider the features of the solar and 
the interplanetary magnetic fields in view of 
above results. (b) Next, we discuss the behavior, 
in the solar atmosphere, of the electron com- 
ponent of solar low energy cosmic rays. And 
then, we obtain some results that the energies 
of the electron component are order of 100Kev 
—l1 Mev, and that these electrons will lose 
their energy by the synchrotron radiation in 
the sunspot magnetic fields (100 Gauss-1000 
Gauss) extending in the solar outer atmosphere, 
and will not propagate into outer space. From 
this view point of these calculation, we are 
able to examine the acceleration mechanism and 
propagation one of the low energy cosmic rays 
produced from the sun. 


SAKURAI K. (Geophysical Institute, Kyoto 


University) The Trapped Orbits in the 
Stormer Theory We consider the trapping 
mechanism of charged particles in the earth’s 
dipole magnetic field in view of Stormer. theory. 
(1) As we can treat the motion of a charged par- 
ticle in the axisymmetric magnetic field as a 
two dimensional problem, we are able to apply 
this method to the case of the dipole magnetic 
field and obtain the general forms of the equa- 
tions of motion of a charged particle. (2) Since 
the forms of this problem are equivalent to the 
motion of a particle in a kind of potential field, 
we are able to examine the general properties 
of motion of this particle without solving the 
equations of motion. (3) In consequence, we 
obtain some results that there are two ways 
with regard to the trapping mechanism of 
charged particles as shown in the following; 
(a) periodic trapped orbits, (b) temporarily 
trapped orbits. It is possible to apply the above 
consideration to explain the cause of Van Allen 
radiation outer belt. Thus, we can examine 
the rigidities of the trapped particles in this 
belt and the stability of this trapping mechanism. 


AONO Y. and Koizumi T. (Radio Re- 
search Laboratories, Japan) Study on Strati- 
fication of F Region For the purpose of 
improving the F1 scaling practices, an investi- 
gation has been made on stratification of the 
F region of the ionosphere, using the daily f- 
plot focussed in 1958 for several stations in 
various latitudes so as to make clear the latitude 
effect. The stations, of which the f-plot is used 
in this study, are indicated in geomagnetic 
latitude as follows: (i) High latitude group 
Ottawa (N 56.9°), Moscow (N 50.8°), White Sands 
(N 41.2°), (ii) Middle latitude group Yamagawa 
(N 20.3°), Okinawa (N 15.5°), and (iii) Low 
latitude group Talara (N 6.6°), Chiclayo (N 5.3°), 
The results of this analysis give the conclusion 
describe below. The stratification of F region 
is so stable in high latitude in summer season, 
that the median and quartile values of the 
frequency at stratification seem meaningful, 
while they will be meaningless in every seasons 
or in any latitudes other than in high latitude 
in summer season because of the unstability of 
stratification. In general, the stratification be- 
comes stable in every latitude, when foF2 at 
noon is less than 8 Mc. According to the in- 
vestigation on the daily f-plot for various sta- 
tions, most of the frequency at stratification 
appears in the value of (5+1.5) Mc from 10h 
to 14h around noon. Therefore it might be 
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called the F1 layer for the stratification, at 
which the frequency changes in this frequency 
band. If the stratification appears sometimes 
in higher or lower frequency band, these strati- 
fication might be called F 1.5 or F 0.5 for ex- 
ample. The diurnal variation of the frequency 
at stratification is composed of such a repeti- 
tion with period of about ten minutes that the 
frequency decreases in a sequential way and 
increases suddenly on appearance of a new 
stratification. Therefore, the geophysical in- 
vestigation on the stratification of F region will 
be made only in use of the f-plot observed in 
every five minutes, at least on quarter-hourly 
basis. The monthly median value of the frequ- 
ency at stratification will be affected by solar 
activity in a similar way as in foF2 but not 
in a simple way of cosine value of the solar 
zenith angle as in fok. In other words, it is 
noticeable that the monthly median value of fo- 
Fl at a certain local day time has a semi- 
annual variation. The value of h’F does not 
change even if the stratification made rapid 
and sequential changes. The value of h’F2 
makes in gradual and regular change of a large 
scale independing on the change of stratification 
of F region. 


MATSUMOTO H. (Institute of Electronics, 
Kyoto University) A Model of the Diurnal 
Variation of NV (A) Profiles in the Ionos- 
phere and its Conductivities A model of 
the diurnal variation of electron density profiles 
with height is constructed from the recent ex- 
perimental data and the theoretical considera- 
tion for the purpose of the study of ionospheric 
physics. References are made to the following ; 
(1) Above the height of # region maximum 
electron density (hm/)—Informations deduced 
from the measurement of Faraday rotation by 
O.K. Garriott. (2) Below the hmF—Tables of 
ionospheric electron density published by the 
Cavendish Laboratory. (3) Above the hmE (E-F 
valley)-Results of rocket measurement in U.S. 
and U.S.S.R. (4) Below the hm£ a) A theoreti- 
cal consideration checked by the results of im- 
pulse sounder and rocket measurements. We 
adopt a modified Chapman layer with the scale 
height with a linear increase vs. height i.e. 
H=Ha-+T(h—ha), and especially, the ion pro- 
duction rate used here in the case of glazing 
incidence is 

I(x,h)=Imnoexp(1+T) {1 —Z—g(y,h)e7}, 
where 


ath F x teh, ath 
FERN) — i sin)’ | Le 


athe siny - (14%) 
lca jolts \] cosec? Ada, 


1=7 logs, Ino=1(0, hmo); 
Ano= atl(hmo—ha); 

a: earth radius. We use a value of [=0.2 
for the E region according to the recent 
results deduced by Kallmann from rocket and 
satellite measurements. b) A variation of the 
night time £& region is based on the theoretical 
works by S. K. Mitra, and checked by the 
results of rocket measurements and _ those 
experimentaly obtained at Kyoto University. 
With this model, variations of electrical 
conductivities o), 6, 02, Gre, xy, Syy versus 
time and height are calculated for the iono- 
sphere in the middle latitude. 


HAKURA Y._ (Hiraiso Radio Wave Obser- 
vatory, Radio Research Laboratories), NAGAI 
M. and SANO K._ (Kakioka Magnetic Obser- 
vatory) On the Relation between Magnetic 
Storm and Ionospheric Disturbances on 
Sept. 13th 1957 An investigation was made 
on the relation between the world-wide patterns 
of ionospheric storms and corresponding current 
systems of geomagnetic storm for each stages 
of a severe upper atmospheric disturbance on 
Sept. 13th, 1957. A close relation of the polar 
magnetic disturbances to the occurrences of polar 
blackouts and storm Es was certified, sugges- 
ting that the current of polar magnetic storm 
flowed on the conductive region of disturbed 
lower ionosphere. It was also found that polar 
cap current system was produced several hours 
before the SC of geomagnetic storm, coinciding 
with the growth of polar cap blackout. After 
the SC, it developed its region to form the 
current system known as that of initial phase. 
At the main phase, ionospheric storm progressed 
to the world-wide scale; the regions of. polar 
blackouts and storm Es elongated towards lower 
latitudes as the Dst field developed. And it was 
shown that auroral currents of bay-like distur- 
bances at that stage, flowed exactly over such 
disturbed regions of high conductivity. 


SINNO K. _ (Hiraiso Radio Wave Observa- 
tory, Radio Research Laboratories) Some 
Characteristics on Solar Particles which 
Excite Abnormal Ionization in the Polar 


Upper Atmosphere A model on the ejec- 
tion and propagation mechanisms which excite 
abnormal ionization in the polar upper atmos- 
phere is presented from some statistical con- 
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siderations of the ionospheric data fmin. The 
evidence that the travelling time (4¢) between 
the major flares and the onset of polar absor- 
ptions show some heliographic relation, is rea- 
sonably interpreted by the present model; that 
is, the particles ejected from the flare are 
trapped into the magnetic cloud which would 
excite geomagnetic storm within 3 days, and 
are carried with rather slow velocity, then, the 
trapped particles could reach to the earth along 
the solar magnetic lines of force which extend 
to the earth, with high velocity. Existence of 
two different onsets of polar abnormal absor- 
ptions, sudden and gradual type of onset, are 
also favorable with this model. 


NaGAlI M., SANO Y. and YANAGIHARA 
K. (Kakioka Magnetic Observatory) On the 
Dst and DS of Each Individual Magnetic 
Storm during IGY (I) For each individual 
storm, the averages of hourly mean value of 5 
low latitude stations widely distributed in longi- 
tudes are used for Dst at equator. At a station 
Dst is considered to be /(@)-Dsto, where Dsto 
is the Dst at equator and /(@) is a function of 
the geomagnetic colatitude 0. For several high 
latitude stations, DS=D—f(@)-Dsto is calculat- 
ed. D’s are hourly values of the disturbance 
field during storm, and f(@) is deduced from 
Chapman-Sugiura’s statistical results. Then 
the ratio A= DS/SDp is calculated for each hour 
of the first 2 days of the storm, where SDo is 
the mean disturbance daily variation by Chap- 
man and Sugiura. Developments and decays of 
Dsto and 4 are considered for each individual 
storm. In this first report, the analyses for 9 
storms in the IGY and some preliminary con- 
sideration of the relation to the other pheno- 
mena, for example, activity of pt pulsation are 
described. 


YAMAGUCHI Y. (Kakioka Magnetic Ob- 
servatory) S.i. Phenomena in the High 
Latitude In the previous paper, this author 
indicated the four foundamental types of s.i., 
based on the horizontal component traces. And 
the impulsive decreasing of the horizontal com- 
ponent occurred as frequently as the increasing 
impulse at Kakioka. In the case of ssc, the 
decreasing of the horizontal component is only 
rare in the middle or low latitude and mainly 
observed in the high latitude. The variations 
in the high latitude of the impulsive decreasing 
in the middle or/and low latitude are examined 
and resulted in that the decreasing components 


do not perfectly depend on the local time. The 
ssc are said to be expressed as follows ; 
Dst+Ds 

and the decreasings of the horizontal component 
in the cases of reversed ssc are observed in the 
region where the horizontal components of Ds 
are negative. Thus, in spite of the resembrance 
of s.i. and ssc, the theory of ssc treating the 
sudden increase of the horizontal component 
does not seem to be suited in the original des- 
cription. From the view point, some theories 
of ssc are reviewed. 


YAMAMOTO M. and MAEDA H. (Ge- 
ophysical Institute, Kyoto University) The 
Simultaneity of Geomagnetic Sudden Im- 
pulses ——The start time of geomagnetic sudden 
impulses was examined by using IGY data. 
Analysis of rapid-run magnetograms at more 
than ten widely-separated magnetic stations 
yielded following results: (1) Time differences 
of sudden impulses around the earth were within 
one minute. (2) The sudden impulses always 
occurred first in high latitudes. (3) The average 
propagation velocity of sudden impulses between 
Honolulu and Koror was about 1300 km/sec. 
These results are very similar to those obtained 
by Williams for sudden commencements of 
magnetic storms. 


YOSHIMATSU T. (Kakioka Magnetic Ob- 
servatory) Universal Time Daily Inequality 
of the Time of Maximum Depression of 
ssc in Storm-Time (II) In the previous 
paper (Memoirs of the Kakioka Magnetic Ob- 
servatory, 1959) the author showed that the 
time interval between sc and the maximum 
depression of the horizontal intensity, Li: 
undergoes a remarkable universal diurnal in- 
equality against the beginning time of sc, ex- 
pressed by the nearest hour in UT. From the 
point of view of this diurnal inequality it was 
also pointed out that some storms may be com- 
posed of two or more storms, and a few storms 
have abnormally shorter T,’s than would be 
expected, and so on. In this paper some other 
characteristics of 7, are further examined and 
reported as follows: (1) T, does not show any 
diurnal inequality aganist both the times of the 
beginning of the main phase and the last 
phase. (2) The amplitude of the diurnal in- 
equality of JT, is of comparable order to the 
mean of 7;, and the ratio of the maximum 
range of 7; to the mean value seems to become 
larger with decreasing geomagnetic latitude. 
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(3) The times of duration of the initial phase 
show the similar universal time diurnal in- 
equality as that of 7,;, though they are more 
scattered about the hourly mean values than 
those of the latter. (4) The universal diurnal 
inequality of T is closely connected to the dis- 
tribution of the geomagnetic non-dipole field, 
especially in lower geomagnetic latitudes at the 


subsolar meridian corresponding to the universal 
time of sc. ; 


MAEDA H. and ONDOH T. 
cal Institute, Kyoto University) Evidence of 
the Ionospheric Propagation of Hydro- 
magnetic Waves Caused by Nuclear Ex- 
plosions over Johnston Island On ex- 
amining carefully the induction magnetograms 
at Onagawa and Shimosato, Japan, we found 
geomagnetic changes which seemed to be due to 
high-altitude nuclear explosions over Johnston 
Island carried out on August 1 and 12, 1958. 
The changes started 10 sec later than the ex- 
plosion time and lasted out about 30 sec. It 
seems from these results that the changes found 
at stations in Japan were due to hydromagnetic 
waves propagated in the upper ionospheric region 
across the lines of geomagnetic force from the 
explosion point. If so, the propagation velocity 
is estimated at about 500 km/sec for the first 
arrival and about 100km/sec for the last arrival, 
on the average of both days. Moreover, if we 
assume that the initial changes at Honolulu 
were also caused by hydromagnetic waves, the 
average damping distance is estimated at about 
1600 km. 


(Geophysi- 


MAEDA H., SAKURAI K. and YAMA- 
MOTO M. (Geophysical Institute, Kyoto Uni- 
versity) Solar Radio Waves and the Earth 
Storms——Recently, the studies on the relation- 
ships between the solar phenomena and the 
earth storms have been accumulated rapidly. 
Especially, it seems that the studies with regard 
to the relation of the solar radio emissions to 
the earth storms (cosmic ray storms, ionospheric 
storms and geomagnetic storms) will give a 
effective clue to the studies of the cause and 
the mechanism of these storms. We examine 
systematically this problem, with regard to 
various data during the periods of IGY and IGC, 
and, referring to the studies done already, try 
to explain the solar-terrestrial relationships 
either unificatively or quantitatively. At first, 
in Part (1), we investigate the following pro- 
blems; (1) the relation of the geomagnetic 
storm producing solar flares and solar radio out- 


bursts type II and IV, (2) the relationships 
among the geomagnetic storms, the associated 
cosmic ray increases, and the solar radio out- 
bursts of type IV. 


NAGATA T., OGUTI T. and THOMATSU 
T. (Geophysical Institute, University of 
Tokyo) Polar Storms in July 1959 
Very close relation between geomagnetic vari- 
ation and simultaneous luminosity of aurora was 
found in case of geomagnetic bay, giant pulsa- 
tion, rapid pulsation and ssc* observed at Syowa 
station (lat. 69°.0S, long. 39°.6E). Horizontal 
disturbance vectors € of the geomagnetic field 
and zenith luminosity J of auroral green line, 


2 5577A, are mainly examined, with reference 
of maximum electron density of sporadic E 
whenever possible. The results show that 
generally 

J=K ®. 
In case of isolated geomagnetic bay, the amount 
of coefficient K is 3x10-° KR/y. In case of 
pulsations, however, periodic variations dé and 
4jJ are superposed on respective back ground &) 
and J), and consequently 

4J=k4é, 
where k is in rough agreement with 2 K&. 
Relation between 4/ and 4& in case of ssc* is 
also in agreement with the above mentioned 
general conclusion. Therefore, this may give a 
direct evidence for that ssc*, or more generally 
polar part of ssc, is caused by impinging cor- 
puscles in the #-region of the ionosphere accom- 
panied by auroral displays. As for rapid pulsa- 
tion of period less than 10 seconds, 4é is under 
effect of impedance of electric conductivity ow- 
ing to recombination of electrons and ions, 
expressed by 

n=Q—a'n?, 

where Q and a’ denote respectively ionization 
rate and effective recombination. Since a’~10~* 
cgs in the auroral zone ionosphere, Q~a’n® 
holds for variations of period of longer than 
several minutes. But for variation of period of 
10 sec or less, the impedance effect becomes 
significant. By considering this effect, 4 //4&~7 
KR/y for &)=2007 in case of rapid variation can 
be well explained in comparison with 4 J/4é= 
1 KR/y for the longer period variations. 
An exceptional case was found for geo- 
magnetic giant pulsation, where A/c = x0 
KR/y for &~0. This pulsation of period of 
about 6 min may be called the second type 
giant pulsation, and denoted by Pg II. In com- 
parison with that geomagnetic bay, giant pulsa- 
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tion of type I, rapid pulsation and ssc* are 
originated in the /-region of the ionosphere 
accompanied by auroral displays of reasonable 
luminosity, Pg II is accompanied by very weak 
aurora. This result would indicate that Pg II 
is originated in outside of the ionosphere. 


MAEDA RR. (Faculty of Science, Tokyo 


University) and TAKEUCHI H. (The Institute 
of Physical and Chemical Research) Temporal 
Variation and Spatial Distribution of 
Radiation at 1000 km above the Earth 
near Japan——Temporal variation and spatial 
distribution at 1000 km above the earth near 
Japan are investigated by analyzing the data 
received at Itabasi, Tokyo from U.S. Satellite 
1959 Iota during March and April in 1960. 
From April 6th, to April 13th great fluctu- 
ations of radiation occurred above the moderate 
latitude region, but were not seen above the 
low latitude. By cosmic ray decrease and geo- 
magnetic storm on March 31 th, radiation inten- 
sity did not vary so much above the lower latitude 
as above the higher latitude. Spatial distribu- 
tion of radiation at 1000 km above the moderate 
latitude is different from that above the higher 
latitude. The former has stationary contours 
of radiation intensity which are nearly parallel 
to the geomagnetic latitude lines, but the latter 
has irregularities and fluctuations. In order to 
interpret the variation and the distribution of 
the radiation, the outer zone and the slot of 
Van Allen Belt are worthy of consideration. 


MATUURA N. (Geophysical Institute, 
Tokyo University) On the Earth Storms (II) 
Interaction between the Solar Corpuscular 
Stream and the Earth’s Magnetic Field — 
The general features of electric current caused by 
charged particles trapped in the earth’s mag- 
tic field have been examined. There are two 
parts in the electric current; one is caused by 
drift motions of trapped charged particles owing 
to the nonuniformity of the magnetic field and 
the other is caused by the particles of which 
guiding centres are lying outside the small 
volume under consideration and of which gyra- 
tion paths are cutting the small volume. There- 
fore, the net electric current is given by the 
sum of the two. The latter part of electric 
current becomes complex when the physical 
quantities, such as density, magnetic field and 
particle velocity, vary sharply within the length 
‘of the order of Larmor radius of the particle. 
This problem has been discussed also. It is 


impossible to interpret the cause of main phase 
of the magnetic storm by a simple consideration 
of drift motion of the trapped charged particles. 
Further, the effect of temporal variation of 
magnetic field on trapped charged particles has 
been examined. 


TOHMATSU T. (Geophysical Institute, 
Tokyo University) On the Earth Storm 
(II) Energy and Flux of the Corpuscular 
Stream Impinging the Earth’s Atmosphere 
(II) The simultaneous ionization and optical 
excitations by energetic electrons in the aurora 
are investigated by considering an isotropically 
incident electron flux whose energy spectrum is 
inversely proportional to the kinetic energy in 
the range between 100 eV and 100 KeV, as is 
consistent with the rocket results. The statio- 
nary electron density profile and the lumi- 
nosities of N,* negative bands and [OI] 5577 
radiation are calculated by using available ioni- 
zation and excitation cross sections. The results 
are compared with those from the optical and 
radio observations. It is concluded that the 
energetic electrons as well as the protons are 
able to be the main ionizing and exciting source 
in aurorae. The appearance of N, positive 
groups and its local time variation are discussed 
in view of the spatial and time variations of 
the energy spectrum of impinging stream. The 
intensity of the positive groups relative to the 
total luminosity becomes a usefull measure of 
the higher energy component of the incoming 
flux. 


OcuTI T. (Geophysical Institute, Tokyo 
University) On the Earth Storm (II) Inter- 
relations among the Upper Atmosphere 
Disturbance Phenomena—— Statistical rela- 
tions among the polar elementary magnetic 
storms including bay disturbances, giant puls- 
ations, and short period pulsations, polar iono- 
spheric disturbances, for example, abrupt in- 
crease in f)9Es and time change of h’Es, and time 
dependence of height of the auroral lower border, 
were examined. The occurrence frequency of 
positive, sharp negative and broad negative bays 
takes the maximum value centered around 18, 23 
and 3h LT respectively, corresponding to the 
three maxima of f,Es and of auroral appearance 
in high latitudes. Most of giant pulsations, how- 
ever, occur in morning-daytime site (6-14h LT) 
of the auroral zone, having no remarkable 
related disturbance in the ionospheric and auro- 
ral phenomena, and may therefore be attribu- 
table to the geomagnetic disturbances propagated 
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from the exosphere. Worldwide patterns of 
h’Es and auroral lower border are quite similar 
to each other. The former which takes the 
maximum height of about 140 km at 18-20h LT 
in the auroral zone, gradually depressed through 
the night, and finally attains to the minimum 
hight of about 100 km at 3h LT. It rises again 
at 6h to about 120km. From the facts men- 
tioned above, it is concluded that the height of 
disturbances represented by anomalous ioniza- 
tion and simultaneous excitation in the lower 
ionosphere in the auroral zone is depressed 
systematically with respect to local time. The 
conclusion is supported also by the systematic 
increase of appearance probability of 1 P.G. in 
the auroral spectra thruogh the night. The 
local time dependence of the disturbance height 
may be due to the change in the sort and/or 
in the primary energy of impinging corpuscles 
during the polar earth storms. 


KOKUBUN S. (Geophysical Institute, Tokyo 
University) On the Earth Storm (II) Geo- 
magnetic Disturbances in the Polar Re- 
gions (I1)——Characteristics of individual geo- 
magnetic disturbances at the geomagnetic con- 
jugate points in the northern and southern 
polar regions are examined based upon the IGY 
data. i) One to one correspondence of isolated 
bay-type variations between the geomagnetic 
conjugate points is especially good for case of local 


night time variations with typical auroral display. 


Correlation coefficient of horizontal disturbing 
force at Baker Lake (73.7°, 315.1°, geomagnetic 
coordinate) and Little America (74.0°, 312.0°, 
geomagnetic coordinate) amounts to 0.85. This 
result may indicate that the corpuscular stream 
originated in the exosphere tends to flow in 
along geomagnetic lines of force toward both 
their north and south ends in the ionospheric 
regions. In case of large magnetic storm, cor- 
relation of disturbing force between these two 
points becomes much poorer. This may suggest 
that relation between two points linked by a 
geomagnetic line of force is much disturbed 
owing to turbulent field in the outer atmos- 
phere. It was found that giant pulsations of 
period 4"~15” occurs simultaneously with 
nearly same period at the geomagnetic con- 
jugate points in the local morning. It may be 
worthwhile to note further that bay-type varia- 
tion in local night time has a tendency to shift 
westward to geomagnetic line of force. The 
velocity of this drift is 70°~150°/hr in order 
of magnitude. 


KANEDA E, (Geophysical Institute, Tokyo 
University) On the Earth Storms (II) 
Photometric Studies of Aurora All-Sky 
Camera Records As a new analytical 
method, ASCA-photoelectric analyser (ASCA- 
P.A.) has been designed for the quantitative 
analysis of aurora all-sky camera records. By 
the photometric method, ASCA-P.A. transforms 
aurora luminosity in a choiced area of movie- 
projected image of aurora all-sky records into 
deflection by pen writing oscillograph. Though 
the stage of analysis by ASCA-P.A. remains a 
preliminary one, future prospect of analysis by 
ASCA-P.A. is very promising. 


YANAGIHARA K,  (Kakioka Magnetic 
Observatory) Observation of Geomagnetic 


Pulsations of Very Short Period, 0.1 to 
10 sec Design of measuring apparatus of 
geomagnetic pulsations having very short period, 
0.1 to 10sec, is outlined. Preliminary observa- 
tions made at Kakioka by the apparatus, show 
some interesting results including the inter- 
mittent enhancements of activity during storms 
or bay-type disturbances. In particular, the 
pulsations of period, 2 to 4 sec, are predominant. 
On the other hand, irregular noises are ob- 
served at Kakioka with appreciable amplitudes 
in the geomagnetically calm period by the high 
sensitive apparatus. Stray currents of artificial 
origin, for example the electric railway near 
Tokyo, bring the noises, maybe. More accurate 
observation of pulsations by the apparatus will 
be made soon after at Memambetsu where arti- 
ficial noises are reduced. 


UTASHIRO S. (Hydrographic Division, 
Maritime Safety Agency) On the Local 
Character of the Geomagnetic Pulsation 
Pe The author studied the results observed 
by induction magnetometers at the Japanese 
four stations (Memanbetsu, Onagawa, Simosato 
and Kanoya). From the results of the studies, 
local inequality of Pe was found. Period or 
mode of Pc change with latitude, and amplitude 
of Pc becomes larger with increasing latitude. 
Such local Pe type pulsation may be produced 
by the hydromagnetic oscillations in the region 
between the inner Van Allen Belt and iono- 
sphere. Two mode of oscillations exist in this 
region, toroidal and poloidal oscillations. In the 
higher latitude the poloidal and toroidal oscil- 
lations are observed, in the lower latitude only 
the poloidal oscillation is observed. Therefore, 
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local Pe occurs frequently in the higher latitude 
than in the lower latitude. Then, the author stud- 
ied on the local character of Pc observed at equa- 
tor and lower latitude during the International 
Geophysical Year by rapid run magnetometer, 
that is, at the Koror, Guam, Apia and Honolulu 
Magnetic Observatory. Comparing the data of 
these observatories, Pc type pulsation occurs 
simultaneously at equator and lower latitude, 
and phase and period is the same at each ob- 
servatory. From these phenomena, it seems 
that only Pc type pulsations produced by the 
poloidal part of the hydromagnetic oscillation 
are observed at equator and lower latitude. 


KAWAMURA M., OuUSHIMA  H., 
Kurusu K. and YANAGIHARA K. (Kaki- 
oka Magnetic Observatory) Wide Distribu- 
tion of the Disturbance Vector of Pulsa- 
tion Pe Disturbance vectors of pc having 
the period of about 20 sec are measured on the 
microfilms of the quick-run magnetograms for 
the 10 stations, Lo, Mb, Kr. Gu, Ho, PB, Co, 
Si, Tu and Fr. Worldwide development and 
decay of disturbance field for 24 hours of April 
29 1958 are considered. Enhancement of am- 
plitude in the day hemisphere is clearly noticed 
as deduced from the statistical daily variation 
of activity at one station. But the vector dis- 
tribution changes so irregularly from time to 
time that it is very difficult to suppose an 
equivalent overhead current system. 


SATO T. (Ist Research & Development 
Center, Japan Defence Agency) A Morphology 
of Geomagnetic Giant Pulsation——Some 
statistical results of characteristics of geomagnet- 
ic giant pulsation that has period of about one 
to ten minutes are shown using the magneto- 
grams during the IGY. Giant pulsation appears 
in the stations over the world, generally between 
3” and 2" L.T.. The amplitude and occurrence 
frequency increase with increasing geomagnetic 
activity, especially, in great magnetic storm 
both are amplified. The results of the analysis 
of some great SC storms show that the giant 
pulsation is a world-widely simultaneous phenom- 
enon, and at a certain universal time the 
periods and phases of the horizontal component 
are in coincidence with each other over the 
world. The results also show that the pulsation 
has a close relation with a bay disturbance in 
northern and southern hemispheres and the 
pulsation on the calm day in one hemisphere 
corresponds to the bay disturbance in other 


hemisphere. From these facts, the giant pulsa- 
tion may be caused by the propagation of the 
hydromagnetic wave originated near the earth 
when the charged particles responsible for the 
bay disturbance impinge into the upper atmos- 
phere. The micropulsation with a period of 10 
to 40 seconds has also close correlation with a 
bay disturbance and the frequency of occur- 
rence, amplitude dependency with latitude are 
similar with those of pulsation mentioned above. 
Therefore the origin of the former seems to be 
same as that of the latter. 


Kato Y. (Geophysical Institute, Tohoku 


University) Geomagnetic Pulsation Accom- 
panying with the Solar Flare Usually 
we cannot detect the micropulsation caused by 
the solar flare in the record of the induction mag- 
netograph, but at the time of very intense solar 
flares, experienced during IGY and also suc- 
ceeded year, we observed very remarkable geo- 
magnetic pulsations, continued in a few minutes 
having the period of 70-80 sec. There is a 
possibility that some kind of hydromagnetic 
oscillation in the earth’s outer atmosphere will 
be excited by a sudden increase of solar radia- 
tion accompanying with the solar flare. If the 
intensity of external radiation, illuminating a 
highly conducting gas which rotates rigidly with 
the magnetic field, increases suddenly and the 
radiative equilibrium condition is established 
during a very short time in a rotating gas, then 
the torsional oscillation of line of force will be 
excited through the magnetic restoring force. 
Therefore it is probable to consider that above 
observed geomagnetic pulsation (especially it 
is very intense on the E-W component) is 
caused by these torsional oscillation of the outer 
atmosphere excited by the solar flare. 


SAITO T. (Geophysical Institute, Tohoku 
University) Geomagnetic Pulsations in the 


North American Region—Damped Type 
Pulsations Accompanying with ssc 
Using the rapid-run magnetograms in the 
North American region during IGY, fairly 
damped type pulsations accompanying 
with ssc are analyzed. These types of pul- 
sations are divided into two groups as follows. 
(1) Pulsations with period of a few hundred 
seconds. At high latitude, ssc often accom- 
panies with a damped type pulsation in the 
period range 200 to 300 sec. The amplitude of 
the pulsation increases with geomagnetic lati- 
tude and reaches its maximum at a little inside 
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of the auroral zone. (2) Pulsation in the period 
range 10 to 40 sec. Damping time of this type 
of pulsation is only a few minutes, nevertheless, 
it is very similar with pe type pulsation in the 
following characters. (i) Period contains is 
longer in the daytime than in the night time. (ii) 
Amplitude is larger in the daytime. (iii) Am- 
plitude is also larger in the higher latitude. 
Then, a suggestion is given that the resonated 
hydromagnetic oscillation in the lower part of 
the exosphere gives rise to the above pulsation 
or pc, especially the former is excited by the 
hydromagnetic shock wave of ssc. 


SINNO K. (Hiraiso Radio Wave Observatory, | 


Radio Research Laboratories) An Analysis 
on the Micropulsations from the Data of 
Wide Spread Observations——Some results 
on the analysis of the micropulsations are pre- 
sented from wide spread data of rapid run 
magnetograms during IGY. Pc micropulsation 
is divided into two frequency bands which is 
shorter than 30 seconds and longer than 40 
seconds, shorter Pe and longer Pc, respectively. 
And, also, existence of another pulsation with 
3 to 7 minutes period which is closely related 
with the longer Pc, is confirmed. Pf micro- 
pulsation can be divided into two groups: viz. 
Pt- which is a Pt associated with a negative 
bay in the auroral zone and P?t* which isa Pt 
either associated with a positive bay or one 
which appears with no bay in the auroral zone. 
The equivalent current systems of the micro- 
pulsation Pc’s and Pi’s, are also derived. In 
this regard, a tentative discussion is given of 
hydromagnetic oscillations in the outer atmos- 
phere are excited by solar corpuscular streams. 


TAMAO T. (Geophysical Institute, Tohoku 


University) Deviation of Geomagnetic Field 
and Hydromagnetic Characteristics in the 
Outer Exosphere With relation to hydro- 
magnetic interactions between charged particles 
and the disturbing geomagnetic field, prelimi- 
nary considerations on the mechanism of cur- 
rent system responsible to the deviation of geo- 
magnetic field from the dipoler one which ob- 
served at 5~7 earth’s radii by Sonnet et al. 
are given. In order to account this deviation, 
gas pressure of the order of 10-° dyne/cm sec 
is required. Pressure expected only from ther- 
mal particles is not sufficient and the extention 
in space of deviation of field resulting from the 
consideration of pressure balance is too large 
to explain the observed results. Betatron 


acceleration of charged particles in the weak 
hydromagnetic disturbing region is discussed 
and the energy spectrum of particles character- 
ized by the threshold velocity of the order of 
10° cm/sec obtained. It may be likely that the 
source of the observed deviation of field is a 
ring current resulting from particles of few 
Kev diffused into and accelerated within the 
region of 7~12 earth’s radii. 


NAMIKAWA T. (Institute of Polytechnic, 
Osaka City University) Propagation of 
Hydromagnetic Wave in the Earth’s Upper 
Atmosphere The laws of reflection and 
refraction of hydromagnetic wave in the earth’s 
upper atmosphere were found, using Pidding- 
ton’s equation of motion of partially ionized 
gases. 


KAMIYAMA H. and OkuDA M. (Geo- 
physical Institute, Tohoku University) Inten- 
sity Variation of [OI] 6300A Emission in 
Night Airglow——It has often been pointed 
out that the intensity of the red oxygen lines 
is closely related to the maximum electron den- 
sity in the / region of the ionosphere. Intro- 
ducing many parameters as well as maximum 
electron density, we carried out more reasona- 
ble calculation. Oxygen red lines in night air- 
glow have been suggested to be produced 
through 

s k 
O*+0,—>0-+0,+ 
a 
O,* +e—>O’/+0” (1) 


O’—>0 + hy 
the rate coefficient of which are denoted. by k 
and a, respectively. Then, the intensity of red 
lines observed on the ground is given by 


i= al 7rn{0.*] -nfeldz, (2) 


where Z, is the lower boundary of the F region 
and Z, is an appropriate altitude. As n [O,*] 


k 
can be approximated to 7102], we can calcu- 


late the above equation. If we assume that the 
distributions of [O,] and x [e] are given by 
n{O,]=n,[0,]-exp{—A(Z—Z,q)}, 
Zi: datum level (3) 
1/0: scale height 
and 
2 
nfel=nntel{1-— "+, 
ym 
ym: semi-thickness of the layer (4) 
(Z>=Zm—Ym) 
we have 
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(5) 
when Z,=Zn+¥n=Z1 
2knm 1 
2 Shatel {nf +m[0:]—7,— (m{0.] 


ny [0.1], (6) 


where 2) and #, are the values at Z and Z,, 
respectively. Examining the observational data 
of the airglow and of the ionosphere, we found 
out the set of the highest correlation coefficients, 
being about 0.8 on the average, between them 


1 
when we abopt the following values : @ ookm 


1 
and Z,=Z, in winter, and en km and Z, 


=300km in summer. We could also show that 
the observed results of the airglow intensity 
were given by Jj,,=/+/, where Jy is nearly 
constant in a course of a night. But J) exhibits 
the marked seasonal variation, and is suggested 
to be contributed from the lower part F region 
excluded from the parabolic ionosphere assumed 
by (4). 


NAKAMURA 1. (Tokyo Astronomical 
Observatory) Latitude Effect of Oxygen 
Red Line of Night Airglow and their 
Relation to Ionospheric F Layer——The 
photoelectric observations of the emission lines 
5577A, 5893A, 6300A and OH bands in the in- 
fra-red in night airglow were made on board 
the m/s Soya in the 4th JARE. The follow- 
ings are the results reduced from the mea- 
surement of 6300A red-line. (1) The intensity 
of the red oxygen line shows remarkable in- 
crease in low latitude and also have high season- 
al variation. (2) In low latitude, daily varia- 
tion of the intensity is larger than in high 
latitude. It seems that this effect is due to the 
vertical motion of the emission layer in low 
latitude. (3) The absolute intensity R of the 
red line shows good correlation with Q which 
is expressed by 


Q= k\* n(e)+n(O,)dz 
1 
where n(e) is the electron density in the Fy 
layer, 7(O,) is the density of oxygen molecules, 
and K is the recombination factor of O, and 
O*. From this fact we may conclude that the 


following reaction predominate in emitting 6300A 
line. 


O++0,—>0,*+0 
O,*-+e—>O*+O0** 
O*(D)—>O(?P2) + hy(6300) 
The value of K becomes to about 1.3-10-!! cm*/ 
sec in the F layer. 


NAKAMURA M. (Institute for Optical 
Research, Tokyo University of Education) 
Measurements of Altitude of the Glow- 
ing Layer of OH Infrared Bands in the 
Night Sky (I) Trigonometrical Method 
——tThe glowing layers in the night sky are, 
in general, neither homogeneous nor stationary, 
but display strong variation spacially and with 
time. If any patches are observed simultaneous- 
ly at two stations situated at an appropriate dis- 
tance, one can find the altitude of the glowing 
layer by means of a trigonometrical method. 
The OH emission in the night airglow was 
observed at Kakioka, Ibaragi Prefecture and 
Maruyama, Chiba Prefecture, which situate 
about 135km apart in almost north and south. 
Simultaneous observations were succeeded in 
three nights up to date. (Dec. 27/28, 28/29, and 
29/30, 1959). Isophoto maps on which the in- 
tensity distribution of OH emission was shown 
by contour lines were drawn for each 30 minutes, 
and about 30 pairs of simultaneous isophoto 
maps were obtained. The altitudes are estimat- 
ed at 74km from the most probable shift of 
the simultaneous maps, 72 km from the intensity 
on the same meridians, and higher than 50 km 
and lower than 100 km from the time variations 
of the mean intensities of superposed area. 


ONAKA R., NAKAMURA M. and TAKE- 
ZAWA S. (Institute for Optical Research, 
Tokyo University of Education) Measurements 
of Altitude of the Glowing Layer of OH 
Infrared Bands in the Night Sky (II) 
Estimation by Van Rhijn Formula with 
Data taken at a High Altitude by Jet 
Plane To reduce the atmospheric extinction 
of light from the night airglow, measurements 
were carried out on a jet plane at a_ high 
altitude (elevation was about 13km). Seven 
zenith distances—0°, 50°, 65°, 70°, 75°, 80° and 
85° were surveyed. Six flights were carried 
out and measurements were performed about 
15 times for each flight. The altitude of the 
glowing layer of OH bands was estimated by 
comparison with Van Rhijn formula at 60 km 
from the earth surface. 


WAKAI N. (Radio Research Laboratories) 
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The Movement of Hs Cloud and Night F 
Layer at Syowa Base The movement of 
Es cloud; Using the sequences of ionograms at 
Syowa Base, the drift velocity of Es cloud was 
deduced from the change of the apparent height 
of oblique Es trace when it was assumed that 
Es cloud was moving horizontally toward or 
away from the station. Observed curves of ap- 
parent height versus time are classified roughly 


into two groups of different velocity of the motion. 


First group has a horizontal velocity of about 
100 meters per second which is identical with 
the value obtained from the other method to 
measure the wind velocity in E region. On the 
other hand, second group has a velocity of 800 
meters per second. It is noticed that the 
amount of the velocity is similar with the 
change of the position to be bombarded by 
ionizing agent. Night FE layer; Several character- 
istics of night E layer were analyzed and com 
pared with those at other high latitude stations. 
Night £& layer could be often observed prior to 
an increase of Es ionization and more frequent- 
ly observed at night time even in winter when 
sun is below horizon at noon. 


NAKAMURA J. (College of General Edu- 
cation and Tokyo Astronomical Observatory, 
University of Tokyo) Spectra of Aurorae (I) 
Time and Space Statistics Aurorae 
may be classified into three types; (a) yellow- 
greenish one, (b) red lower border and (c) red 
upper border or red all over. Representative 
spectral lines of them are as follows; (a) OLA 
5577 line and N, N.G., (b) Ist P.G. and Ha and 
(c) OJ2A 6300, 6363 lines. Space statistics were 
performed as to aurorae which were appeared 
at Showa Base (69°00S, 39°35’E) during 1959. 
Auroral maximal zone lies in (a) southern region 
of our base, (b) northern region of there and (c) 
somewhat southern bearings of there respective- 
ly. Their height may be as follows; (a) about 
100 km, (b) 80~90 km and (c) abut 200~300 km. 
Time statistics: (a) type aurora appears at daily 
night time and has nearly constant local time 
sequences. (b) occurs with magnetic bay which 
maximum activity occurs near midnight. Always 
Hea line appears first and one hour or two later 
Ist P.G. appears. (c) occurs accompanied by 
magnetic storms and it seems to have no cor- 
relation with local time. 


NAKAMURA J. (College of General Edu- 
cation and Tokyo Astronomical Observatory, 
University of Tokyo) Spectra of Aurorae 


(II) Doppler Shift of Ha Line—Dispersion 
of auroral spectra taken at Showa Base was 
320A per mm and it could be read to about 
1.5A by Hartmann’s formula. All of Ha and 
Hf lines were suffered Doppler shift, and quan- 
tity of them varied by each direction along 
magnetic meridian. If we assumed the incident 
direction is coinside with the direction of ma- 
gnetic lines of force (azimuth 350°, zenith angle 
26°), above variation is well elucidated. The 
approaching velocity along lines of force dis- 
tributes between 1.0~9.5 x 10? cm per sec, which 
maximum lies about at 3.0x107cm per sec. 
A velocity component perpendicular to approa- 
ching one is about one third of the approaching 
velocity that is distributed between 0.5~3.0x 
10’cm per sec. 


ONDOH T., KITAMURA T. and MAE- 
DA H. (Geophysical Institute, Kyoto Univer- 
sity) The Relation of the Dawn Chorus 
to the Ionospheric Absorption of Cosmic 
Noise and the Auroral and Geomagnetic 
Activities We have studied the relation 
of the dawn chorus to the ionospheric absorp- 
tion and the auroral and geomagnetic activities 
on geomagnetically quiet and disturbed days 
from September, 1959 to March, 1960, by using 
the data observed at College, Alaska. There isa 
negative correlation between the auroral activity 
and the chorus. This fact seems to imply a possi- 
bility that the chorus is absorbed by unusual 
ionization due to auroral particles. On the 
other hand, the correlation between the chorus 
and the cosmic noise is positive. This may 
mean that both the chorus and the cosmic noise 
are influenced by ionospheric absorption in a 
similar manner. If the origin of the chorus 
and the cosmic noise is independent of magnet- 
ic storms, there may be expected a functional 
relation between them. The origin of the chorus 
seems, however, to be affected by magnetic 
storms, so that the deviations from such a 
relation may be regarded as an effect of magnet- 
ic storms. These deviations are in a negative 
correlation to the magnetic K-index. This fact — 
may imply that the exospheric region where the 
geomagnetic line of force starting from the 
earth’s surface at 64° in geomagnetic latitude 
goes up on magnetically quiet days becomes 
inactive in generation of the dawn chorus with 
the increase of magnetic activity. Such a ex- 
planation will be clarified by using more data 
at other stations. 
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ONDOH T. (Geophysical Institute, Kyoto 
University) On the Origin of the V.L.F. 
Noise in the Earth’s Exosphere—tThe ra- 
diation power of the Cerenkov radiation and the 
proton cyclotron radiation are respectively 1.7 
x 10-4 watts. m~ and 7.2x10-* watts. m~? at 
2 earth-radii from the center of the earth on the 
geomagnetic line of force passing through the 
earth’s surface at the geomagnetic latitude 60°, 
where the electron density is 5x10*/cm® and 
the velocity of the protons is 10®°cm/sec. As 
Ellis has pointed out, the Cerenkov radiation 
due to charged particles moving along the 
geomagnetic line of force gives the same type 
frequency characteristics as that of the travell- 
ing wave tube type amplification given by 
Gallet. The condition of the generation of the 
Cerenkov radiation and the operation of the 
travelling wave tube type amplification is given 
by V, cosd=— and Viz respectively, where 
V, and V, are the velocities of protons. When 
@ is small but not zero, we may consider that 
the whistler mode radio waves radiated by the 
mechanism of the Cerenkov radiation are am- 
plified by the amplification mechanism of the 
travelling wave tube type due to different charg- 
ed particles at the different part in the space 
and the time. At 2 earth-radii from the center 
of the earth on the geomagnetic line of force 
starting from the earth’s surface at the geo- 
magnetic latitude 60°, we have for the total 
emission in a column of length 200 km which 
consists of fast protons (density 10™4/cem*),,. 1.4 
x10-* watts. m-® (c/s)"!._ The received power 
flux of dawn chorus signals is of the order of 
10-* to 10-'* watts. m~? (c/s)-!._ Thus we can 
not obtain theoretically the received power flux 
without the gain of amplification greater than 
60 db. 


INOUE Y. and HAsHIzuME S. (Group 


of Radio-physics, Electrical Engineering Divi- 
sion, Meguro Research Laboratories, National 
Defence Agency) Exospheric Radio Sound- 
ing——Two new methods studying the exo- 
sphere are presented. (1) Whistler Compres- 
sion Filter (WCF) ; A whistler is generated by a 
lightning impulse. Its dispersion is formed in 
the course of propagating along the dispersive 
delaying path in the exosphere. Its delaying 


D(Gm) 


time is fg= er for short whistler, where 


D(6@m) is the dispersion at the magnetic latitude 
Om. After a whistler receiver, a phase 
equalizer with the frequency response of 
exp {27i(2D(0)/./ 7 )f} is followed: This equal- 
izer compensates the dispersive delaying time 
of a whistler and converts it into a sharp pulse 
like an original lightning inpulse. WCF does 
not enhance noises and disperses sferics. It 
has a high discrimination for a whistler. By 
WCF, the 24 hr clock-run observation becomes 
possible. (2) Inverse Long Whistler (ILW); 
A signal is transmitted with the inverse fre- 
quency dispersion to a long whistler. The 
signal is converted into a sharp pulse by 
the dispersive delay in the course of propaga- 
tion in the exosphere. The signal equivalent to 
a lightning impulse in the greatest class is as 
follows: 

transmitting power.............. 3KW 

transmitting duration............ 1 sec 

frequency dispersion range .... 1-10KC 
The echo-pulse can be detected with an ample 
S/N. Furthermore, by adjusting the frequency 
dispersion of the signal, we can know precisely 
the characters of the propagation of VLF-radio 
waves in the exosphere. ILW can be trans- 
mitted in sferics-calm times. By ILW, it is 
possible to study the propagation exactly along 
the geomagnetic line of force. As its PRF is 
some 0.25c/s in middle latitudes, one can in- 
vestigate the rapid variations of the exospheric 
characteristics. 
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